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Executive  Summary 


The  U.S.  Army  Corps  of  Engineers  software  HEC-FDA  is  a  tool  that  uses  Risk-based 
Analysis  to  compute  Expected  Annual  Damage  (EAD)  for  flood  damage  reduction 
studies.  EAD  reduction  is  computed  as  the  difference  between  EAD  with  and  without 
alternative  projects,  a  quantity  used  to  aid  in  flood  damage  reduction  project  selection. 
Depending  on  the  time  and  money  spent  in  the  collection  and  preparation  of  the  input 
data  sets,  relatively  lesser  or  greater  accuracy  can  be  achieved  in  EAD  computation. 
However,  although  the  estimate  can  be  improved,  some  parameter  and  model 
uncertainty  is  unavoidable.  It  is  important  to  examine  the  sensitivity  of  EAD  and 
EAD  reduction  to  realistic  differences  in  input  functions  and  their  associated 
uncertainties.  These  results  can  be  used  to  identify  factors  with  the  least  and  greatest 
impact  on  the  resulting  EAD  to  make  recommendations  on  the  best  investments  to 
improve  estimates  of  EAD  and  EAD  reduction. 

The  current  study  addresses  the  following  objectives  using  HEC-FDA: 

>  Examine  and  compare  the  sensitivity  of  Expected  Annual  Damage  (EAD)  and 
EAD  reduction  to  variations  in  the  three  primary  input  functions,  namely  flow- 
exceedance  probability,  stage-flow  and  damage-stage. 

>  Examine  and  compare  the  sensitivity  of  EAD  and  EAD  reduction  to  uncertainty  in 
the  three  primary  input  functions. 

>  Identify  the  factors  with  the  least  and  greatest  uncertainty  contribution  to  EAD. 
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>  Use  the  results  of  the  analyses  to  recommend  the  best  investments  in  methods  and 
data  to  improve  flood  damage  analysis. 

To  examine  possible  differences  in  sensitivity  and  uncertainty  between  streams  with 
different  hydrologic,  hydraulic  and  economic  characteristics,  34  HEC-FDA  data  sets 
from  across  the  continental  United  States  and  Hawaii  were  sorted  according  to 
drainage  area  and  slope.  Four  categories  were  created  to  form  a  two  by  two  matrix  of 
drainage  area  and  stream  slope.  Based  on  the  information  contained  in  the  data  sets 
according  to  this  matrix,  representative  synthetic  functions  were  created  as  base  cases 
for  a  typical  large  drainage  basin  with  flat  stream  slope  and  a  typical  small  basin  with 
steep  stream  slope. 

For  each  category,  the  base  case  mean  flow  was  taken  from  the  middle  of  the  sorted 
data  sets,  its  standard  deviation  based  on  an  assumed  coefficient  of  variation  of  0.5, 
and  its  skew  assumed  to  be  zero  for  an  assumed  Log  Pearson  Type  III  (LPIII) 
distribution  of  flood  peaks.  An  appropriate  stage-flow  curve  was  generated  based  on  a 
correlation  developed  between  mean  flow  and  stage-flow  curve  average  slope.  A 
damage-stage  curve  was  generated  based  on  typical  curves  in  the  sorted  data  sets  and 
curves  used  in  previous  studies  (Amell,  1989;  Beard,  1990).  The  proportion  of  total 
damage  is  used  rather  than  total  damage  to  put  all  analyses  on  the  same  economic 
basis. 

For  the  sensitivity  analysis,  parameters  were  varied  in  each  of  the  primary  input 
functions  to  calculate  EAD  and  EAD  Reduction  with  50,  100  and  250-year  (0.02,  0.01 
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and  0.004  exceedance  probability)  levees.  Mean,  standard  deviation  and  skew  of  the 
flow-exceedance  probability  function,  y-intercept  and  vertical  location  of  the  stage- 
flow  function,  and  lower  and  upper  bounds  and  inflection  point  locations  of  the 
damage-stage  function  were  varied.  To  examine  sensitivity  to  uncertainty,  equivalent 
record  length  and  standard  deviation  of  error  about  the  stage-flow  and  damage-stage 
functions  were  varied.  To  complete  the  uncertainty  analysis,  the  difference  between 
the  0.25  and  0.75  exceedance  values  from  the  EAD  distribution  were  extracted  from 
the  output  to  measure  the  uncertainty  contribution  from  each  input  function.  Based  on 
the  EAD  calculated  in  each  of  the  model  runs,  elasticity  was  calculated  as  the  percent 
difference  in  EAD  produced  per  percent  change  in  the  parameter  being  varied. 

To  test  the  results  of  the  numerical  experiments  on  the  synthetic  data  sets,  sensitivity 
to  uncertainty  was  examined  for  two  real  sample  data  sets.  A  data  set  from  a  typical 
large  drainage  basin  on  the  Blue  River,  Missouri  and  a  data  set  from  a  typical  small 
drainage  basin  on  the  Chippewa  River,  Georgia  were  used  for  comparison.  The  sample 
data  sets  served  to  verify  results  obtained  from  the  synthetic  data  sets  and  expand  the 
results  to  include  damage-stage  functions  with  different  shapes. 

The  numerical  experiments  showed  significant  differences  between  the  elasticity  of 
EAD  and  EAD  reduction  for  the  large  and  small  basins.  Figure  A  summarizes  the 
elasticity  of  EAD  to  parameters  of  its  major  input  functions,  and  Figure  B  summarizes 
the  elasticity  of  EAD  to  uncertainty. 
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Figure  A.  Sensitivity  Analysis  Summary 

Where  Mean  Flow  =  Mean  flow  of  the  flow-exceedance  probability  function;  S.D. 
Flow  =  Standard  deviation  of  flow-exceedance  probability  function;  Skew  =  Skew  of 
the  flow-exceedance  probability  function;  Stage-Flow  B  =  Average  slope  of  the  stage- 
flow  function;  Stage-Flow  C  =  Y-intercept  of  the  stage- flow  function;  Dam-Stage 
LB/UB  =  Lower  Bound/Upper  Bound  of  damage-producing  stages  in  the  damage- 
stage  function;  Dam- Stage  A/B  =  First  and  second  inflection  points  in  the  damage- 
stage  function. 
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Figure  B.  Sensitivity  to  Uncertainty  Analysis  Summary 


The  results  show  that  EAD  and  EAD  reduction  can  be  relatively  sensitive  to  different 
variables  over  different  ranges.  EAD  reduction  is  more  sensitive  to  changes  in 
function  values  for  larger  levees  than  for  smaller  levees.  In  general,  EAD  and  EAD 
reduction  are  most  sensitive  to  stage-flow  function  slope  (depending  on  the  basin  area, 
shape  of  the  stream  channel  and  floodplain)  and  the  upper  and  lower  bounds  of  the 
damage-stage  function  (depending  on  the  location,  height  and  types  of  structures  in  the 
floodplain),  and  least  sensitive  to  flow-exceedance  probability  skew  (generally  a 
regional  characteristic).  EAD  and  EAD  reduction  are  generally  more  sensitive  to 
uncertainty  in  flow-exceedance  probability  due  to  a  limited  flow  record  length  than  to 
stage-flow  uncertainty  due  to  rating  curve  scatter  or  uncertainty  in  Manning’s  n, 
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except  for  situations  in  a  small  basin  with  shorter  flow  record  lengths  and  larger  stage- 
flow  uncertainty.  It  is  most  important  to  get  an  accurate  estimate  of  the  variables  with 
the  highest  elasticity,  such  as  stage-flow  function  shape,  and  least  important  to  get  an 
accurate  estimate  of  the  variables  with  the  lowest  elasticity,  such  as  flow-exceedance 
probability  skew.  The  value  of  additional  information  can  be  used  in  investment 
decision-making  to  estimate  the  benefits  expected  from  improving  median  input 
function  specification  and  reducing  uncertainty. 
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Chapter  1 
Introduction 

The  U.S.  Army  Corps  of  Engineers  software  HEC-FDA  is  a  tool  that  uses  Risk-based 
Analysis  to  compute  Expected  Annual  Damage  (EAD)  for  flood  damage  reduction 
studies.  EAD  reduction  is  computed  as  the  difference  between  EAD  with  and  without 
alternative  projects,  a  quantity  used  to  aid  in  flood  damage  reduction  project  selection. 
Depending  on  the  time  and  money  spent  in  the  collection  and  preparation  of  the  input 
data  sets,  relatively  lesser  or  greater  accuracy  can  be  achieved  in  EAD  computation. 
However,  although  the  estimate  can  be  improved,  some  parameter  and  model 
uncertainty  is  unavoidable.  It  is  important  to  examine  the  sensitivity  of  EAD  and 
EAD  reduction  to  realistic  variations  in  input  functions  and  their  associated 
uncertainties.  These  results  can  be  used  to  identify  factors  with  the  least  and  greatest 
impact  on  the  resulting  EAD  to  make  recommendations  on  the  best  investments  for 
improving  estimates  of  EAD  and  EAD  reduction. 

1.1  Objectives 

The  current  study  addresses  the  following  objectives  using  HEC-FDA: 

>  Examine  and  compare  the  sensitivity  of  Expected  Annual  Damages  (EAD)  and 
EAD  reduction  to  changes  in  the  three  primary  input  functions,  namely  flow- 
exceedance  probability,  stage-flow  and  damage-stage.  This  will  be  a  sensitivity 
analysis  as  defined  in  (Morgan  and  Henrion,  1990,  pg.  172),  “methods  for 
computing  the  effects  of  changes  in  inputs  on  model  predictions”. 
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>  Examine  and  compare  the  sensitivity  of  EAD  and  EAD  reduction  to  uncertainty  in 
the  three  primary  input  functions. 

>  Identify  the  factors  with  the  least  and  greatest  contribution  to  EAD  uncertainty. 
This  will  be  an  uncertainty  analysis  as  defined  in  (Morgan  and  Henrion,  1990,  pg. 
1 72),  “methods  for  comparing  the  importance  of  the  input  uncertainties  in  terms  of 
their  relative  contributions  to  uncertainty  in  the  outputs”. 

>  Use  the  results  of  the  analyses  to  recommend  the  best  investments  in  methods  and 
data  to  improve  flood  damage  analysis. 

1 .2  Study  Overview 

This  study  examines  the  sensitivity  of  HEC-FDA  results  to  the  following  parameters 
using  a  set  of  numerical  experiments.  The  first  three  sets  of  experiments  vary  the 
shape  of  the  input  functions  analytically,  while  the  second  three  vary  the  specified 
uncertainty.  Under  each  analysis,  the  following  variables  were  varied  directly: 

1 .  Sensitivity  to  Flow  Exceedance  Probability  Function 

Mean  Flow 

Standard  Deviation  of  Flow 
Skew 

2.  Sensitivity  to  Stage-Flow  Function 

Slope  of  Stage-Flow  Function 
-  Y-intercept  of  Stage-Flow  Function 

3.  Sensitivity  to  Damage-Stage  Function 
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Lower  bound  of  damage-producing  stages 
Upper  bound  of  damage-producing  stages 
Location  of  lower  inflection  point  in  function 
Location  of  upper  inflection  point  in  function 

4.  Sensitivity  to  Flow  Exceedance  Probability  Uncertainty  and  Uncertainty  Analysis 
-  Equivalent  Record  Length 

5 .  Sensitivity  to  Stage-Flow  Uncertainty  and  Uncertainty  Analysis 

Standard  Deviation  of  Error 

6.  Sensitivity  to  Damage-Stage  Uncertainty  and  Uncertainty  Analysis 

Standard  Deviation  of  Error 

The  numerical  experiments  listed  above  were  conducted  on  synthetic  base  case 
functions  typical  of  two  general  types  of  damage  reaches  -  a  flat  stream  with  large 
drainage  area  and  a  steep  stream  with  small  drainage  area.  To  test  the  results  of  the 
synthetic  functions,  the  sensitivity  to  uncertainty  experiments  were  conducted  on  two 
actual  sample  data  sets.  Sensitivity  indices  were  developed  to  provide  a  consistent 
basis  for  comparison.  The  results  of  these  numerical  experiments  were  used  to 
address  the  stated  objectives. 
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Chapter  2 
Background 


2.1  EAD  Computation  by  Direct  Integration 

Expected  Annual  Damage  (EAD)  has  been  calculated  traditionally  without  explicit 
incorporation  of  uncertainty.  To  compute  EAD,  a  flow-exceedance  probability 
function  is  generated  from  a  historical  flow  record,  rainfall-runoff  calculations  or 
regional  information;  a  stage-flow  relationship  is  approximated  from  field 
measurements,  water  surface  profile  analysis  or  hydraulic  routing;  and  a  damage-stage 
relationship  is  approximated  from  records  and  surveys.  For  each  exceedance 
probability,  the  corresponding  flow  is  obtained  from  the  flow-exceedance  probability 
curve.  The  median  stage  associated  with  that  flow  is  obtained  from  the  stage-flow 
curve  and  the  median  damage  corresponding  to  that  stage  is  obtained  from  the 
damage-stage  curve.  After  going  through  this  process  for  a  range  of  specified  flow- 
exceedance  probabilities,  a  damage-probability  curve  is  generated.  With  direct 
integration,  the  area  under  the  damage-probability  curve  yields  expected  annual 
damage.  The  difference  in  EAD  between  the  with-  and  without-project  conditions  has 
been  used  to  estimate  the  benefits  expected  from  a  project  to  aid  in  project  selection. 
This  process  is  illustrated  in  Figure  1  (US  ACE,  1989). 


10 


FLOW- EXCEEDANCE  PROBABILITY 


STAGE-FLOW 


Figure  1.  EAD  Computation  with  Direct  Integration 


2.2  Expected  Probability  Estimator 

The  first  efforts  to  include  uncertainty  in  flood  frequency  and  flood  damage  analysis 
resulted  in  the  development  of  the  expected  probability  estimator  (Beard,  1960,  1978; 
Stedinger,  1983b).  The  expected  probability  estimator  was  developed  to  produce  a 
better  estimate  of  the  true  probability  of  occurrences  over  a  broader  scale  (i.e., 
national)  than  a  median  estimator.  For  each  quantile,  the  expected  probability 
estimator  provides  an  estimate  of  the  mean  probability  of  exceeding  a  flow  level.  The 
earliest  work  in  developing  the  expected  probability  estimator  was  by  Beard  (1960), 
and  refined  in  Guidelines  For  Determining  Flood  Flow  Frequency,  Bulletin  #17B  of 
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the  Hydrology  Subcommittee.  The  most  current  revised  methods  for  computing  an 
expected  probability  for  hydrologic  data  are  detailed  in  (IACWD,  1982). 


2.3  Bayesian  Inference 

Bayesian  inference  has  been  advocated  as  a  framework  to  incorporate  uncertainty  in 
federal  investment  decisions  and  flood  frequency  analysis  for  almost  thirty  years  (Al- 
Futaisi  and  Stedinger,  1999;  Bodo,  1976;  Davis,  1972;  Freer  and  Beven,  1996; 
Krzysztofowicz,  1983;  Kuczera,  1999;  Stedinger,  1983;  Wood,  1975).  Bayesian 
inference  is  a  method  of  examining  risk  by  considering  a  set  of  possible  future 
statistics  rather  than  only  the  statistics  of  the  limited  sample.  In  general,  Bayes 
Theorem  estimates  the  posterior  distribution  based  on  a  prior  distribution  and 
likelihood  function.  Mathematically,  Bayes  Theorem  can  be  summarized  as  the 
following: 

f(0|  x)  =  f(0)  '  l(x|  0)/  l(f(0)  '  l(x|  0)  d0) 

Where  l(x|  0)  is  the  likelihood  of  observation  sequence  x,  f(0)  is  the  prior  distribution, 
and  f(0|  x)  is  the  posterior  distribution  (Davis,  1972). 

Limitations  of  Bayesian  inference  for  forecasting  include  the  informativeness  of  the 
prior  distribution  and  the  computational  effort  required.  As  informativeness  of  the 
prior  distribution  decreases,  Bayesian  analysis  can  become  suboptimal  unless 
procedures  are  used  to  account  for  an  uninformative  prior  distribution.  In  general,  the 
informativeness  of  the  prior  distribution  depends  on  whether  or  not  it  is  a 
representative  sample  of  the  population  of  values  (Krzysztofowicz,  1983). 
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The  likelihood  function  generally  represents  the  analyst’s  degree  of  belief  that  a 
parameter  has  different  values  (Al-Futaisi  and  Stedinger,  1999).  Using  the  prior 
distribution  and  likelihood  function,  the  parameters  of  the  posterior  distribution  are 
usually  found  by  simulation.  Posterior  distribution  parameters  can  be  obtained  by 
sampling  from  the  likelihood  function  about  the  prior  distribution.  During  simulation, 
the  posterior  distribution  should  become  more  peaked  at  the  population  value  as  the 
number  of  samples  becomes  large.  Computational  methods  can  be  set  up  to  allow 
updating  of  the  prior  distribution  as  information  is  obtained  (Freer  and  Beven,  1996). 

Krzysztofowicz  (1983)  demonstrates  with  Bayes  Theorem  that  decision-makers  that 
ignore  uncertainty  in  probabilistic  forecasts  always  incur  an  opportunity  loss  as  the 
number  of  observations  increases.  Given  that  the  amount  of  flow  data  and  the  number 
of  federal  investments  in  flood  control  are  expected  to  continue  to  increase,  the  value 
of  applying  Bayesian  Inference  to  flood  damage  analysis  should  increase.  In  light  of 
this  need  to  incorporate  uncertainty  in  flood  damage  analysis  in  a  systematic  fashion, 
Risk-based  Analysis  has  been  under  development. 

2.4  Risk-Based  Analysis  Computation 

In  the  evaluation  of  flood  control  projects,  the  variety  of  uncertainties  makes  it 
difficult  to  determine  system  performance  under  floods  of  varying  magnitude.  The 
USACE  guidelines  give  the  following  operational  definition  of  uncertainty  (NRC, 
1995): 
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Uncertainty:  Uncertain  situations  are  those  in  which  the  probability  of 
potential  outcomes  and  their  results  cannot  be  described  by  objectively  known 
probability  distributions,  or  the  outcomes  themselves,  or  the  results  of  those 
outcomes  are  indeterminate. 

For  flood  damage  analysis,  uncertainty  cannot  be  described  perfectly  by  objectively 
known  probability  distributions.  However,  uncertainty  can  be  included  explicitly  in 
the  analysis  by  assuming  a  best-fit  probability  density  distribution  to  describe  the 
range  of  likely  functions.  The  parameters  of  the  assumed  distribution  can  be 
approximated  from  knowledge  of  the  system,  inferences  from  regional  similar  gaged 
basins  and  data  collection.  The  U.S.  Army  Corps  of  Engineers’  Risk-based  Analysis 
incorporates  uncertainty  with  best-fit  distributions,  as  detailed  below. 

2.4.1  COMPUTATIONAL  METHOD 

The  U.S.  Army  Corps  of  Engineers  has  incorporated  uncertainty  into  flood  damage 
analysis  with  the  computer  software  HEC-FDA.  HEC-FDA  uses  a  Monte  Carlo 
simulation  to  sample  the  interaction  among  the  hydrologic,  hydraulic  and  economic 
relationships  and  their  individual  uncertainties  (USACE,  1996).  The  primary 
functions  are  flow-exceedance  probability,  stage-flow  and  damage-stage.  Each 
relationship  has  a  distribution  of  uncertainty  about  its  median  function,  estimated  from 
knowledge  of  the  system  and/or  data  collection  and  analysis  methods.  Random 
sampling  from  the  range  of  likely  flow-exceedance  probability  distributions,  with  the 
corresponding  range  of  likely  stage-flow  and  damage-stage  probability  distributions  is 
repeated  until  the  confidence  limits  of  annual  damages  are  not  changed  significantly 
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when  more  samples  are  taken.  This  sampling  process  results  in  an  expected  value  of 
annual  damage  (HEC,  1998).  A  graphical  summary  of  this  method  can  be  seen  in 
Figure  2.  The  probability-density  distribution  of  uncertainty  about  each  function  is 
shown  with  the  input  functions  on  the  left,  a  sample  iteration  is  shown  in  the  middle 
and  the  resulting  sample  damage-probability  function  on  the  right.  The  area  under  the 
sample  damage-probability  curve  produces  one  EAD  sample  for  each  iteration.  The 
average  of  all  samples  produces  a  mean,  or  expected  value  of  annual  damage. 


Figure  2.  HEC-FDA  Monte  Carlo  Analysis 
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2.4.2  FLOW-EXCEEDANCE  PROBABILITY  FUNCTION 


The  Water  Resources  Council  published  “Bulletin  17B,  Guidelines  for  Determining 
Flood  Flow  Frequency”  to  promote  a  consistent  approach  to  flood  frequency  analysis. 
In  this  publication,  the  Council  recommends  the  use  of  the  Log  Pearson  Type  III 
(LPIII)  distribution  for  flow  frequency  (IAWCD,  1982).  Lindquist  (1995)  examined 
the  effect  of  choosing  different  distributions  to  describe  the  flow  frequency 
relationship  and  showed  that  the  choice  of  distribution  can  have  a  significant  effect  on 
EAD.  The  LPIII  describes  the  flow-exceedance  probability  distribution  with  three 
parameters,  estimated  using  the  method  of  moments.  By  adjusting  these  three 
moments,  the  LPIII  is  flexible  enough  to  fit  a  wide  variety  of  historical  data  sets 
(Lindquist,  1995).  In  HEC-FDA,  data  sets  that  are  not  described  well  by  LPIII  can  be 
specified  graphically  (HEC,  1998).  The  graphical  methods  are  commonly  used  for 
regulated  flow  conditions.  However,  for  the  purposes  of  setting  up  functions  that  can 
be  altered  systematically  for  sensitivity  analysis,  the  LPIII  distribution  was  used  for 
the  present  study. 

2.4.3  HYDROLOGIC  UNCERTAINTY 

Hydrologic  uncertainty  is  represented  by  the  distribution  of  error  about  the  flow- 
exceedance  probability  function.  The  main  source  of  hydrologic  uncertainty  is  the 
limited  historical  flow  data  available  at  gaging  stations.  Most  gaged  locations  have 
from  10  to  70  years  of  data,  not  all  of  which  are  completely  reliable.  Often,  gages  are 
washed  out  during  the  highest  flow  events,  further  limiting  the  accuracy  of  flood  flow 
readings.  In  addition,  there  is  greater  uncertainty  associated  with  the  measurement  of 
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high  flow,  low  frequency  events  and  estimation  of  their  associated  exceedance 
probabilities  than  with  more  frequent  flows  (HEC,  1998).  For  an  LP1II  distribution  of 
flood  peaks,  hydrologic  uncertainty  is  included  in  flood  damage  analysis  using  a  large 
sample  approximation  to  the  non-central  t  error  distribution  depending  on  gage  record 
length  (IAWCD,  1982;  HEC,  1998).  The  USACE  gives  guidelines  for  approximating 
the  equivalent  record  length  for  gaged  and  non-gaged  locations.  These  guidelines  are 
summarized  in  Table  1  (USACE,  1996b). 

Table  1.  Equivalent  Record  Length  Guidelines 


Method  of  Freauencv  Function  Estimation 

Equivalent  Record  Length1 

Analytical  distribution  fitted  with  long-period  gauged  record  available  at  site 

Systematic  record  length 

Estimated  from  analytical  distribution  fitted  for  long-period  gauge  on  the  same 
stream,  with  upstream  drainage  area  within  20%  of  that  point  of  interest 

90%  to  100%  of  record  length  of 
gauged  location 

Estimated  from  analytical  distribution  fitted  for  long-period  gauge  within  same 
watershed 

50%  to  90%  of  record  length 

Estimated  with  regional  discharge-probability  function  parameters 

Average  length  of  record  used  in 
regional  study 

Estimated  with  rainfall-runoff-routing  model  calibrated  to  several  events 
recorded  at  short-interval  event  gauge  in  watershed 

20  to  30  years 

Estimated  with  rainfall-runoff-routing  model  with  regional  model  parameters 
(no  rainfall-runoff-routing  model  calibration) 

10  to  30  years 

Estimated  with  rainfall-runoff-routing  model  with  handbook  or  textbook  model 
Darameters 

10  to  15  years 

1  Based  on  judgement  to  account  for  the  quality  of  any  data  used  in  the  analysis,  for  the  degree  of  confidence  in 
models,  and  for  previous  experience  with  similar  studies. 

(USACE,  1996b) 


For  normal  and  log  normal  distributions,  the  non-central  t  distribution  has  been 
recommended  for  its  flexibility  in  describing  the  sampling  error  distribution 
conditionally  upon  position  along  the  flow-exceedance  probability  curve  (IAWCD, 
1982;  Stedinger,  1983c;  Tung,  1987;  Chowdhury  and  Stedinger,  1991;  Afshar,  et  al., 
1994;  Lindquist,  1995).  A  study  by  Bao,  et  al.  (1987)  concurs  with  these 
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recommendations  and  offers  an  extension  to  better  incorporate  skew  in  the  uncertainty 
estimate. 

The  effects  of  hydrologic  uncertainty  on  EAD  were  studied  by  Bao,  et  al.  (1987)  for 
normal,  log  normal,  Pearson  Type  III  (PHI),  and  LPIII  flow-exceedance  probability 
distributions.  The  results  of  this  study  show  that  EAD  is  very  sensitive  to  hydrologic 
uncertainty,  particularly  for  record  lengths  of  less  than  60  years.  Results  from 
Stedinger  (1997)  showed  that  damage  estimates  are  significantly  more  sensitive  to 
hydrologic  uncertainty  for  record  lengths  of  less  than  50  years. 

Although  uncertainty  due  to  a  limited  flow  record  is  explicitly  incorporated  in  HEC- 
FDA,  uncertainty  due  to  other  factors  is  more  difficult  to  quantity.  Some  studies  have 
been  completed  to  assess  the  relevance  of  long  term  climate  change  to  flood  control 
project  planning.  Venkatesh,  et  al.  (1999)  found  that  climate  risk  is  of  approximately 
the  same  importance  as  uncertainty  in  the  damage-stage  relationship  in  terms  of  the 
penalty  suffered  if  it  is  ignored  for  Lake  Erie,  a  lake  with  a  large  drainage  area. 
However,  Mendelsohn  (1997)  found  that  because  basin- specific  changes  in  runoff 
from  global  warming  are  currently  uncertain  and  much  delayed,  most  project  analyses 
will  be  unaffected  by  global  warming.  Additional  uncertainty  is  introduced  if  flood 
control  dams  change  operations  in  the  future  or  have  surcharge  capacity  that  was  not 
included  in  the  flood  routing  calculations.  The  amount  of  time  flood  water  remains  in 
the  flood  plain  and  recession  characteristics  are  also  uncertain  due  to  variable  drainage 
characteristics  and  human  efforts  (NRC,  1995).  Although  a  method  has  not  yet  been 


devised  to  explicitly  incorporate  all  of  these  uncertainties,  they  are  important  to 
understand  potential  differences  between  model  predictions  and  reality. 

2.4.4  STAGE-FLOW  FUNCTION 

The  stage-flow  function  is  defined  graphically  in  HEC-FDA  based  on  field 
measurements  and/or  hydraulic  calculations.  Although  there  can  be  significant  scatter 
in  field  measurements,  the  function  is  generally  defined  based  on  a  single  line  of 
increasing  stage  with  flow  (Westphal,  et  al.,  1999).  However,  this  approximation 
becomes  rougher  with  increasing  flow,  as  is  discussed  below. 

2.4.5  HYDRAULIC  UNCERTAINTY 

Hydraulic  uncertainty  is  derived  from  the  estimated  relationship  between  flow  and 
stage.  Turbulence  characteristics  create  variations  in  stage  for  a  given  flow  while 
eroding  the  stream  channel  and  changing  the  cross  section  (NRC,  1995).  In  light  of 
this,  higher  flows  possess  greater  hydraulic  variability  (and  thus  uncertainty)  than 
medium  and  lower  flows.  Westphal,  et  al.  (1999)  demonstrate  that  the  uncertainty 
introduced  by  unsteady,  high  flows  was  significant  for  the  1993  flood  along  the 
Mississippi  River  at  St.  Louis,  Missouri.  These  results  also  imply  a  tendency  towards 
hysteresis  for  extreme  conditions,  or  stage-discharge  relationships  varying  between  the 
rising  and  falling  hydrograph  limb.  However,  the  results  of  the  Mississippi  River 
study  also  found  that  neither  changes  in  stage  nor  changes  in  flow  rate  correlate  with 
short-term  scour  and  fill  at  gauging  sites  (Westphal,  et  al.,  1999). 
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Additional  uncertainty  in  the  stage-flow  relationship  is  introduced  by  variation  in 
flow-stage  measurements,  ice,  debris,  bulking,  seasonal  variations  in  Manning’s  n 
value,  structural  integrity,  multiple  possible  failure  modes,  and  surveying  inaccuracies 
of  levees  and  water  control  structures  (Burnham,  pers.  com.,  2000).  Settling  of  levees 
and  subsidence  of  flood  plain  lands  over  time  also  can  be  a  factor.  Simplifications 
made  for  hydraulic  calculations  such  as  averaging  flow  resistance  characteristics  and 
limiting  the  number  of  cross  sections  and  their  associated  survey  measurement  errors 
also  contribute  to  hydraulic  uncertainty  (NRC,  1995).  Guidelines  for  estimating  a 
standard  deviation  of  error  of  the  hydraulic  uncertainty  distribution  depending  on 
survey  methods  are  given  in  HEC  (1986)  as  part  of  a  detailed  analysis  on  the  sources 
of  stage-flow  error.  The  study  found  that  uncertainty  in  Manning’s  n  is  the  greatest 
contributor  to  stage-flow  error,  although  it  is  not  apparent  exactly  how  it  affects  the 
variance  of  that  error  (HEC,  1986). 

The  best-fit  distribution  of  stage-flow  uncertainty  can  vary  with  how  the  relationship 
was  estimated.  The  scatter  about  the  best-fit  stage-flow  curve  estimated  at  a  gaged 
reach  can  be  approximately  normal  or  more  skewed.  Freeman  et  al.  (1996)  found  that 
the  gamma  distribution  can  represent  a  wide  range  of  rating  curve  error  distributions 
from  normal  to  highly  skewed,  and  recommended  it  to  describe  stage  uncertainty 
(USACE,  1996b).  Methods  for  approximating  the  uncertainty  for  ungaged  stream 
reaches  based  on  measurable  stream  parameters  and  computed  water  surface  profiles 
are  detailed  in  USACE  (1996b).  For  ease  of  sensitivity  analysis  with  HEC-FDA,  the 
normal  distribution  was  used  to  approximate  stage-flow  uncertainty,  with  the  range  of 


20 


standard  deviations  of  error  determined  based  on  Manning’s  n  reliability  (USACE, 
1996b). 

2.4.6  DAMAGE-STAGE  FUNCTION 

The  damage-stage  function  is  calculated  as  a  graphical  function  based  on  procedures 
described  in  the  “National  Economic  Development  Procedures  Manual  -  Urban  Flood 
Damage”,  March  1988,  “National  Economic  Development  Procedure  Manual  -  Urban 
Flood  Damage  -  Volume  II:  Primer  on  Surveying  Flood  Damage  for  Residential 
Structures  and  Contents”,  October  1991,  “Catalog  of  Residential  Depth-Damage 
Functions”,  May  1992,  and  “Analysis  of  Non-Residential  Content  Value  and  Depth- 
Damage  Data  for  Flood  Damage  Reduction  Studies”,  April  1996  (HEC,  1998).  The 
general  shape  of  the  damage-stage  function  usually  varies  between  a  straight  line  and 
S-shaped  curve,  depending  on  land  use  characteristics  (Amell,  1989).  The  most 
common  shape  in  the  U.S.  is  most  similar  to  the  “logistic”,  or  “S”  shape  from  Amell 
(1989),  with  the  steepest  part  of  the  curve  between  the  100-year  and  1,000-year  stages 
(Beard,  1990).  Al-Futaisi  and  Stedinger  (1999)  used  a  quadratic  equation  to  describe 
the  damages  above  a  certain  threshold.  A  quadratic  equation  was  not  used  for  this 
study  due  to  the  difficulty  in  systematically  applying  changes  in  function  shape  to 
represent  a  realistic  range  of  conditions  to  be  examined  for  the  sensitivity  analysis. 

2.4.7  DAMAGE  UNCERTAINTY 

Uncertainty  associated  with  the  damage-stage  relationship  can  be  calculated  based  on 
the  individual  uncertainties  of  its  major  components.  Currently,  HEC-FDA  can 
account  for  uncertainty  in  structure  value,  content  to  structure  value  ratio,  depth- 
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percent  damage  and  surveyed  elevations.  These  uncertainties  are  combined  into  an 
estimated  normal,  triangular,  or  log  normal  distribution  of  error  about  the  mean  (HEC, 
1998).  The  standard  deviation  of  this  error  distribution  is  usually  specified  as  a 
percentage  of  the  damages  (Moser,  1993).  However,  Lindquist  (1995)  showed  that  the 
distribution  of  damage  error  does  not  affect  the  calculation  of  EAD  as  long  as  the 
mean  function  is  specified.  The  distribution  of  uncertainty  is  only  significant  if  the 
distribution  of  EAD  is  under  examination. 

The  error  distribution  does  not  account  for  any  changes  in  velocity  (and  water  surface 
elevation)  as  the  floodwater  spreads  onto  the  flood  plain,  or  the  accumulation  of  debris 
that  can  affect  water  stage  and  damage  at  different  locations.  The  distribution  also 
does  not  account  for  any  differences  in  water  stage  between  the  surrounding  flood 
plain  and  individual  structures.  The  way  in  which  the  duration  of  flooding  affects 
damage  is  also  not  considered  (NRC,  1995).  However,  a  study  by  Afshar,  et  al. 
(1998)  assumed  that  the  effects  of  flooding  time  and  volume  are  minor  for 
structure/building  flood  damage  estimation.  Additional  factors  that  influence  flood 
damage  are  described  in  USACE  (1996b). 

Although  these  sources  of  uncertainty  are  difficult  to  quantify  for  inclusion  in  the 
distribution  of  uncertainty,  they  contribute  to  differences  between  the  model 
approximation  and  reality.  Although  it  might  seem  ideal  to  try  to  quantify  and  include 
all  conceivable  sources  of  uncertainty  in  its  distribution,  it  is  important  to  realize  that 
increasing  complexity  does  not  necessarily  yield  an  increase  in  accuracy  (Pilgrim, 
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1986).  Therefore,  it  is  desirable  to  consider  a  best  estimate  of  the  range  of  possible 
distributions  without  breaking  the  distribution  down  into  too  many  small  pieces. 

2.5  Comparison  of  Methods 

While  all  methods  share  the  goal  of  estimating  expected  annual  damage  to  compare 
flood  damage  reduction  alternatives,  Direct  Integration,  Expected  Probability  and 
Risk-based  Analysis  differ  in  assumptions  about  the  most  accurate  way  to  perform 
flood  damage  analysis  under  imperfectly  known  conditions.  More  detailed  discussion 
on  the  controversial  appropriate  use  and  role  of  these  methods  has  been  published  in 
an  ongoing  series  of  papers  (Beard,  1960,  1978,  1990,  1997,  1998;  Thomas,  1976; 
NRCTF,  1978;  Stedinger,  1983;  Amell,  1989;  Gunasekara  and  Cunnane,  1991; 
Rasmussen  and  Rosbjerg,  1991;  Stedinger,  1997;  Goldman,  1997).  While  Direct 
Integration  computes  the  median  annual  damage  without  explicitly  considering 
uncertainty,  Risk-based  Analysis  computes  the  mean  annual  damage  (expected  value) 
with  probability  density  distributions  of  likely  input  functions.  Since  the  expected 
probability  method  estimates  a  mean  damage,  the  method  is  inherent  in  the  risk-based 
approach.  The  higher  value  of  the  mean  relative  to  the  median  in  a  positively  skewed 
uncertainty  distribution  yields  a  generally  higher  EAD.  A  schematic  of  the  positively 
skewed  hydrologic  uncertainty  distribution  is  given  in  Figure  3. 
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FLOW-EXCEEDANCE  PROBABILITY 


Figure  3.  Skewed  Hydrologic  Uncertainty  Distribution  Schematic 


For  Risk-based  Analysis,  the  benefits  of  explicitly  incorporating  uncertainty  in  the 
EAD  calculation  can  be  estimated  using  the  concept  of  the  expected  value  of  including 
uncertainty  (EVIU).  EVIU  is  the  difference  between  the  expected  value  of  the  optimal 
decision  with  uncertainty  considered  and  the  expected  value  of  the  optimal  decision 
with  uncertainty  ignored  (Morgan  and  Henrion,  1990).  For  flood  damage  analysis, 
this  can  be  defined  as: 

EVIU  =  [EADu.wo  -  EADu,w]  -  [EADW0  -  EADW] 
where 

EADu,wo=  Expected  Annual  Damage  including  uncertainty,  without  project 
EADU>W=  Expected  Annual  Damage  including  uncertainty,  with  project 
EADW0  =  Expected  Annual  Damage  without  considering  uncertainty,  without  project 
EADW  =  Expected  Annual  Damage  without  considering  uncertainty,  with  project. 
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Risk-based  Analysis  also  allows  estimation  of  the  value  of  additional  information,  the 
difference  between  the  confidence  in  the  flood  damage  estimate  with  current 
information,  and  the  confidence  in  the  flood  damage  estimate  if  one  more  unit  of 
information  was  obtained  in  any  of  the  primary  input  functions  (Morgan  and  Henrion, 
1990).  An  extension  of  the  value  of  additional  information  is  the  difference  between 
the  current  net  benefits  expected  from  a  project  and  the  net  benefits  expected  if  one 
more  unit  of  information  was  obtained.  The  concepts  of  EVIU  and  value  of  additional 
information  can  be  used  along  with  the  results  of  the  current  study  in  investment 
decision-making.  The  following  section  presents  a  detailed  description  of  the 
approach  used  to  address  the  current  study  objectives. 
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Chapter  3 

Analysis  Approach 


3.1  Base  Case  Functions 

To  examine  possible  differences  in  sensitivity  and  uncertainty  between  streams  with 
different  hydrologic,  hydraulic  and  economic  characteristics,  34  HEC-FDA  data  sets 
from  across  the  continental  United  States  and  Hawaii  were  sorted  according  to 
drainage  area  and  slope.  Four  categories  were  created  to  form  a  two  by  two  matrix  of 
drainage  area  and  stream  slope.  The  distribution  of  data  sets  within  the  matrix  can  be 
seen  in  Table  2. 


Table  2.  HEC-FDA  Data  Set  Matrix 


Stream  Type 

Small  Drainage  Area 
(<100mi2) 

Large  Drainage  Area 
(>100  mi2) 

Flat  Stream  Slope 
(<3  ft/mi.) 

Coulee  des  Cannes,  LA 

Chester  Creek,  PA 

Lafayette/Vermilion  R.,  LA 

Ocmulgee  River,  GA 

Blue  River,  MO 

Chattahoochee  River,  GA 

Flint  River,  GA 

White  Oak  Bayou,  TX 

Willamette  River,  OR 

DesPlaines  River,  IL 

Cumberland  River,  OH 

Tygart  River,  WV 

Sheyenne  River,  ND 

Etowah,  GA 

Tuolumne  River,  CA 

Missouri  River,  MO 

Steep  Stream 

Slope 
(>3  ft/mi. 

Harmon  Canal,  GA 

Wailupe,  HI 

Perry  Creek,  IA 

Strong  Ranch  Slough,  CA 

Deer  Creek,  MN 

Anacostia  River,  VA 

NW  Anacostia,  VA 

Lower  Mission  Cr,  CA 

Antelope  Cr,  NE 

Beargrass  Cr,  KY 

Rio  De  Flag,  A Z 

Warm  Springs,  CA 

Pecan  Bayou,  TX 

Tres  Rios/Gila  River,  A Z 

Murrieta  Creek,  CA 
Scranton/Lackawanna  R.,  PA 

Greens  Bayou,  LA 

White  River,  IN 
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Small  and  large  drainage  areas  were  divided  at  100  square  miles,  and  steep  and  flat 
stream  slopes  were  divided  at  3  feet/mile.  Since  most  of  the  data  sets  fell  into  one  of 
two  categories  (large  drainage  area/flat  slope  or  small  drainage  area/steep  slope), 
representative  synthetic  functions  were  created  as  base  cases  for  these  two  categories. 
Examination  of  cases  in  the  other  two  categories  can  be  included  in  further  study.  The 
following  procedures  were  used  to  generate  representative  base  case  functions. 

3.1.1  FLOW-EXCEEDANCE  PROBABILITY  FUNCTION 

Within  each  category,  the  flow-exceedance  probability  function  for  each  data  set  was 
plotted.  The  range  of  mean  flows  typical  to  the  two  categories  of  interest  can  be  seen 
in  Tables  3a  and  3b. 


Table  3a.  Flow-Exceedance  Probability  Function  of  Known  Large/Flat 

Data  Sets 


Large/Flat  Data  Set 

Illinois  River,  IL 

430 

Lafayette/Vermilion,  LA 

952 

Sheyenne  River,  ND 

1,200 

DesPlaines  R,  IL 

2,555 

Tuolumne  R,  CA 

5,000 

Etowah  R,  GA 

8,283 

Blue  River,  MO 

10,170 

White  Oak  Bayou,  TX 

12,271 

Chattahoochee  R,  GA 

21,702 

Ocmulgee  River,  GA 

26,448 

Tygart  River,  WV 

26,696 

Flint  River,  GA 

31,623 

Cumberland  River,  OH 

96,694 

Willamette  River,  OR 

136,994 

Missouri  River,  MO 

214,042 
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Table  3b.  Flow-Exceedance  Probability  Function  of  Known  Small/Steep 

Data  Sets 


Small/Steep  Data  Set 

Mean  Flow  (cfs) 

Rio  De  Flag,  A Z 

19 

Warm  Springs,  CA 

70 

Beargrass  Cr,  KY 

401 

Lower  Mission  Cr,  CA 

Deer  Creek,  MN 

534 

Harmon/Chippewa,  GA 

707 

Strong  Slough,  CA 

788 

Wailupe,  HI 

949 

Antelope  Cr,  NE 

1,900 

Perry  Creek,  LA 

2,180 

NW  Anacostia  R,  VA 

4,410 

Anacostia  R,  VA 

5,154 

For  consistency  with  Bulletin  17B  recommendations  and  ease  of  sensitivity  analysis, 
base  case  functions  were  created  as  Log  Pearson  Type  III  with  zero  skew,  or  Log 
Normal  (IAWCD,  1982).  The  mean  flow  for  each  base  case  was  chosen  from  the 
middle  of  the  known  data  sets  in  its  respective  category.  For  the  large  drainage 
area/flat  stream,  a  mean  flow  of  25,000  cfs  was  used.  For  the  small  drainage 
area/steep  stream,  a  mean  flow  of  700  cfs  was  used.  The  coefficient  of  variation 
(COV)  is  the  ratio  of  the  standard  deviation  to  the  mean  of  a  distribution.  Linsley 
(1986)  reports  a  typical  range  of  0.3  -  0.8  for  the  COV  of  floods  in  the  United  States 
in  real  space,  while  A1  Futaisi  and  Stedinger  (1999)  and  Landwehr  (1978)  report  a 
typical  range  of  0.5  -  1.0.  Assuming  the  range  stated  by  A1  Futaisi  and  Stedinger 
(1999)  and  Landwehr  (1978),  the  COV  is  between  0.05  and  0.13  in  log  space.  In 
general,  arid  regions  tend  to  have  a  higher  COV,  whereas  wetter  regions  tend  to  have  a 
lower  COV  (Landwehr,  1978).  A  COV  of  0.05  was  assumed  for  each  base  case  in  log 
space  to  represent  a  climate  with  less  flow  variability.  A  COV  of  0.05  yields  a 
standard  deviation  of  0.22  for  the  large/flat  case  and  0.14  for  the  small/steep  case.  The 
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resulting  base  case  functions  are  shown  in  Figure  4.  To  minimize  uncertainty,  the 
equivalent  record  length  was  assumed  to  be  200  years. 


BASE  CASE  FLOW-EXCEEDANCE  PROBABILITY  FUNCTIONS 


Figure  4.  Base  Case  Flow-Exceedance  Probability  Functions 

3.1.2  STAGE-FLOW  FUNCTION 

To  create  base  case  stage-flow  functions  consistent  with  the  flow-exceedance 
probability  functions,  the  stage-flow  functions  from  the  known  data  sets  were  plotted. 
A  quadratic  function  was  fitted  to  the  data  points  for  flows  greater  than  zero.  The 
function  is  connected  via  a  straight  line  to  (0,0).  The  average  slopes  of  each  function 
were  plotted  against  the  mean  flow  from  the  corresponding  flow-exceedance 
probability  function  to  approximate  a  correlation. 

Based  on  the  plot  scatter  for  the  small/steep  case,  the  corresponding  rating  curve  can 
have  an  approximate  slope  of  0.002  to  0.006.  Based  on  the  plot  scatter  for  the 
large/flat  case,  the  corresponding  rating  curve  can  have  an  approximate  slope  of 
0.0001  to  0.0005.  The  correlations  can  be  seen  in  Figures  5a  and  5b.  Based  on  the 
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range  of  y-intercept  values  for  the  appropriate  case  functions,  the  y-intercept  (“C” 
value)  can  vary  between  0  and  5  for  the  small/steep  case  and  between  5  and  1 0  for  the 
large/flat  case.  The  range  of  flows  covered  by  the  stage-flow  function  must  meet  or 
exceed  the  range  of  flows  in  the  flow  exceedance  probability  function.  To  minimize 
uncertainty  in  the  base  cases,  the  standard  deviation  of  error  was  set  to  0.01  feet.  The 
resulting  base  case  stage-flow  functions  are  (See  Figure  6): 

Small/Steep: 

Y  =  -6E-7X2  +  0.004X  +  1 .55 
Large/Flat: 

Y  =  -5E-10X2  +  0.0002X  +  7.89 
where  Y  =  Stage  (ft)  and  X  =  Flow  (cfs). 

For  flows  higher  than  the  apex  of  the  curve,  a  linear  approximation  was  used  with  a 
slope  equal  to  the  slope  of  the  500  cfs  prior  to  the  apex. 
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Mean  Flow  Versus  Stage-Flow  Curve  Avg  Slope 
Steep  Stream,  Small  Drainage  Area 


Mean  Flow  (cfs) 


Figure  5a.  Mean  Flow  -  Stage-Flow  Curve  Average  Slope  Correlation: 
Steep  Stream,  Small  Drainage  Area 


Figure  5b.  Mean  Flow  -  Stage-Flow  Curve  Average  Slope  Correlation: 
Flat  Stream,  Large  Drainage  Area 
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Base  Case  Stage-Flow  Functions 
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Figure  6.  Base  Case  Stage-Flow  Functions 


3.1.3  DAMAGE-STAGE  FUNCTION 

Like  the  first  two  functions,  the  damage-stage  functions  were  created  by  first  plotting 
the  damage-stage  functions  of  the  known  data  sets.  Based  on  the  typical  shapes  of  the 
known  functions  and  studies  by  Amell  (1989)  and  Beard  (1990),  a  piecewise  linear 
function  was  created  to  allow  the  functions  to  be  more  easily  and  systematically 
adjusted  during  the  set  of  runs  necessary  for  sensitivity  analysis.  The  piecewise  linear 
shape  is  most  closely  related  to  the  “logistic”  shape  in  Amell  (1989),  and 
recommended  by  Beard  (1990)  as  the  most  common  curve  shape  for  floodplains  in  the 
U.S.  Assuming  damage  begins  at  the  mean  stage  associated  with  the  2-year  event 
(exceedance  probability  =  0.5)  and  maximizes  at  the  mean  stage  associated  with  the 
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500-year  event  (exceedance  probability  =  0.002),  the  following  typical  function  was 
derived: 

Between  X=0  and  X=LB,  Y=0 
At  X=1/3(UB-LB)  +  LB,  Y=0. 1 
At  X=2/3(UB-LB)  +  LB,  Y=0.9 
At  X=UB,  Y=1 

where  Y  =  Proportion  of  Total  Damage  (between  0  and  1);  X  =  Stage  (ft) 

LB  =  Lower  Bound  Stage  of  Damage  Function;  UB  =  Upper  Bound  Stage  of  Damage 
Function. 


The  proportion  of  total  damage  is  used  rather  than  total  damage  to  put  all  analyses  on 
the  same  economic  basis.  To  minimize  uncertainty  for  the  base  case,  the  standard 
deviation  of  error  was  set  to  0.01  feet.  The  base  case  damage-stage  functions  can  be 
seen  in  Figure  7. 


Base  Case  Graphical  Damage-Stage  Functions 


Figure  7.  Base  Case  Damage-Stage  Functions 
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3.2  Sensitivity  Analysis 

3.2.1  FLOW-EXCEEDANCE  PROBABILITY  FUNCTION 

To  examine  the  effects  of  changing  the  shape  of  the  most  likely  flow-exceedance 
probability  function  on  EAD  and  EAD  reduction,  the  three  parameters  of  the  Log 
Pearson  Type  III  function  were  varied  independently.  Table  4  summarizes  the 
variation  introduced  to  the  flow-exceedance  probability  function  to  examine  the 
sensitivity  of  EAD  and  EAD  reduction  to  flow-exceedance  probability  function 
parameters. 

Table  4.  Flow-Exceedance  Probability  Function  Sensitivity  Analysis 


Base  Case 

Variable 

Base 

Value 

Lower 

Limit 

Upper 

Limit 

Increment 

Large/Flat 

Mean 

25,000  cfs 

12,500  cfs 

25,000  cfs 

500  cfs 

Small/Steep 

Mean 

700  cfs 

350  cfs 

700  cfs 

50  cfs 

Large/Flat 

ismsi 

0.1 

Small/Steep 

EEIHH 

bhh 

mam 

0.37 

0.1 

Large/Flat 

Skew  (log) 

EllMiW 

-0.7 

0.4 

0.1 

Small/Steep 

Skew  (log) 

0 

-0.7 

0.4 

0.1 

The  range  in  mean  values  was  chosen  to  keep  a  reasonable  coefficient  of  variation 


according  to  Landwehr  (1978)  without  having  to  change  the  standard  deviation. 
Accordingly,  the  range  in  standard  deviation  values  was  chosen  to  yield  a  reasonable 
range  in  coefficient  of  variation.  The  range  of  skew  values  was  based  on  the  USGS 
map  of  Generalized  Skew  Coefficients  of  Logarithms  of  Annual  Maximum 
Streamflow  by  One  Degree  Quadrangles  (IACWD,  1982). 


3.2.2  STAGE-FLOW  FUNCTION 

The  base  case  stage-flow  function  was  created  as  a  quadratic  equation  so  its 
coefficients  could  be  varied  systematically.  The  function  is  of  the  form: 
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S  —  aQ2  +  bQ  +  c 

where  S  =  Stage  (ft)  and  Q  =  Flow  (cfs). 

To  vary  the  function  slope,  “b”  was  varied.  To  vary  the  function  vertically,  “c”  was 
varied.  Table  5  shows  the  method  used  to  examine  the  sensitivity  of  EAD  and  EAD 
reduction  to  variation  in  the  stage-flow  function. 


Table  5.  Stage-Flow  Function  Sensitivity  Analysis 


Base  Case 

Variable 

Base 

Value 

Lower  Limit 

Upper  Limit 

Increment 

Large/Flat 

b/Slope 

IMIMI 

0.0001 

0.0001 

Small/Steep 

b/Slope 

gilifiM 

0.002 

0.001 

Large/Flat 

c/Y-Int. 

7.89 

5 

10 

1 

Small/Steep 

c/Y -Int. 

1.55 

0 

5 

1 

It  is  important  to  note  that  the  function  type  assumed  does  not  account  for  reverse  flow 
conditions.  Such  conditions  can  result  from  a  combination  of  backwater  effects  (from 
tides  or  larger  order  streams),  intense  precipitation  and  high  flows.  Reverse  flow 
conditions  are  generally  rare  except  in  flat  streams.  For  example,  this  situation  occurs 
along  the  Vermilion  River  in  Louisiana  (US ACE,  1995).  In  this  case,  it  is  important 
to  have  the  rating  curve  reflect  a  probability  weighted  average  condition  useful  for 
transforming  the  flow-exceedance  probability  curve  to  a  stage-exceedance  probability 
curve  (Goldman,  pers.  comm.,  2000). 


3.2.3  DAMAGE-STAGE  FUNCTION 

The  base  case  damage-stage  function  was  created  as  a  piecewise  linear  function  to 
allow  systematic  variation.  Fitting  polynomial  functions  to  data  points  as  was  done 
for  the  stage-flow  function  proved  inadequate  to  properly  describe  the  damage-stage 
function  shapes  while  coefficients  were  varied.  For  the  piecewise  linear  function, 
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inflection  points  in  the  function  were  specified  according  to  the  following  general 
function: 

Between  X=0  and  X=LB,  Y=0 
At  X  =  A(UB-LB)  +  LB,  Y=0. 1 
At  X  =  B(UB-LB)  +  LB,  Y=0.9 
At  X  =  UB,  Y=1 

where  Y  =  Proportion  of  Total  Damage  (between  0  and  1) 

X  =  Stage  (ft) 

A  =  Location  of  first  inflection  point  (as  fraction  of  total  range  of  damage-producing 
stages) 

B  =  Location  of  the  second  inflection  point 
LB  =  Lower  Bound  Stage  of  Damage  Function 
UB  =  Upper  Bound  Stage  of  Damage  Function. 

For  the  base  case,  A=l/3  and  B=2/3  for  the  flat/large  and  steep/small  cases.  For 
sensitivity  analysis,  the  lower  bound  stage,  upper  bound  stage,  A  and  B  were  varied  to 
examine  the  sensitivity  to  overall  slope  of  the  function  and  slopes  of  three  main 
sections  of  the  function.  Table  6  summarizes  the  method  used  to  examine  the 
sensitivity  of  EAD  and  EAD  reduction  to  differences  in  the  damage-stage  function. 
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Table  6.  Damage-Stage  Function  Sensitivity  Analysis 


Base  Case 

Variable 

Base  Value 

Lower 

Limit 

Upper 

Limit 

Increment 

Large/Flat 

LB 

12.6 

8.6 

16.6 

1 

UB 

23.6 

19.6 

27.6 

1 

A 

1/3 

1/6 

3/6 

1/6 

B 

2/3 

3/6 

5/6 

1/6 

Small/Steep 

LB 

4.2 

2.2 

4.2 

1 

UB 

8.2 

6.2 

8.2 

1 

A 

1/3 

1/6 

3/6 

1/6 

B 

2/3 

3/6 

5/6 

1/6 

LB  =  Lower  bound  of  damage-producing  stages,  UB  =  Upper  bound  of  damage- 


producing  stages,  A  =  Location  of  first  inflection  point  (as  fraction  of  total  range  of 
damage-producing  stages),  B  =  Location  of  second  inflection  point 


3.3  Sensitivity  to  Uncertainty 

In  addition  to  examining  the  sensitivity  of  EAD  and  EAD  reduction  to  differences  in 
the  shape  of  the  input  functions,  its  sensitivity  to  differences  in  uncertainty  was 
examined.  In  HEC-FDA,  uncertainty  can  be  specified  in  each  of  the  three  primary 
input  functions.  Flow-exceedance  probability  uncertainty  is  estimated  as  a  function  of 
the  sampling  error  in  the  mean  and  standard  deviation.  The  uncertainty  in  the  mean  is 
known  to  be  described  by  a  normal  distribution  and  the  standard  deviation  by  a  chi- 
squared  distribution  (Goldman,  pers.  comm.,  2000).  Stage-flow  uncertainty  and 
damage-stage  uncertainty  are  assumed  to  be  described  by  a  normal  distribution.  To 
examine  sensitivity  to  parameter  uncertainty,  the  base  case  functions  were  used  with 
minimal  uncertainty.  Uncertainty  was  estimated  and  introduced  in  increments  for 
each  function. 
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3.3.1  FLOW-EXCEEDANCE  PROBABILITY  UNCERTAINTY 


To  examine  the  sensitivity  of  EAD  and  EAD  reduction  to  flow-exceedance  probability 
uncertainty,  the  equivalent  record  length  was  varied  between  10  years  and  200  years 
for  each  base  case.  As  equivalent  record  length  increased,  uncertainty  and  the 
standard  deviation  of  error  decreased.  For  a  given  equivalent  record  length,  the 
amount  of  uncertainty  increases  with  flow.  Therefore,  the  flat/large  case  possessed 
greater  uncertainty  than  the  steep/small  case  for  the  same  equivalent  record  length. 

As  is  discussed  in  the  previous  section,  the  sampling  error  in  estimates  of  the  mean 
and  standard  deviation  is  used  to  describe  parameter  uncertainty  (Stedinger,  1983b). 
When  skew  is  included  in  the  flow-exceedance  probability  distribution,  the  calculation 
extension  provided  by  Bao,  et  al.  (1987)  is  needed  to  provide  a  better  estimate  of 
parameter  uncertainty. 

Sources  of  error  other  than  a  limited  historical  record  length  can  be  considered  in 
estimating  uncertainty  in  the  flow-exceedance  probability  distribution.  Kuczera 
(1996)  presents  the  results  of  rating  curve  error  on  flood  frequency  inference.  If  the 
flood  flow-exceedance  probability  distribution  is  back  calculated  from  stage 
measurements  and  a  rating  curve,  the  flow-exceedance  probability  distribution 
contains  some  error  in  addition  to  the  error  introduced  by  a  limited  record  length.  The 
error  introduced  by  flood  frequency  inference  can  be  incorporated  into  estimation  of 
equivalent  record  length  as  described  by  USACE  (1996b).  Other  methods  of 
estimating  flow  frequency  such  as  hydrologic  modeling  or  estimation  from  other  gages 
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within  the  watershed  also  can  be  incorporated  into  the  equivalent  record  length 
(USACE,  1996). 

3.3.2  STAGE-FLOW  UNCERTAINTY 

To  examine  the  sensitivity  of  EAD  and  EAD  reduction  to  stage-flow  uncertainty,  the 
standard  deviation  of  error  about  the  mean  function  was  varied.  The  standard 
deviation  of  error  was  assumed  to  increase  linearly  with  flow,  with  zero  uncertainty  at 
(0,0).  The  error  was  assumed  to  remain  constant  above  the  0.01  exceedance 
probability  (100-year)  event.  EAD  and  EAD  reduction  were  calculated  for  standard 
deviation  of  error  equal  to  0.3,  0.6,  1.5  and  3.0  feet  to  represent  varying  accuracy  in 
data.  These  standard  deviations  account  for  likely  ranges  of  error  introduced  by 
survey  technology  limitations,  selected  accuracy,  Manning’s  roughness  coefficient, 
and  stream  hydraulic  properties  for  steady  flow  (HEC,  1986). 

The  effects  of  unsteady  flow  characteristics  become  more  important  at  higher  flows. 
Westphal,  et  al.  (1999)  demonstrate  that  the  uncertainty  introduced  by  unsteady,  high 
flows  was  significant  for  the  1993  flood  along  the  Mississippi  River  at  St.  Louis, 
Missouri.  These  results  also  imply  a  tendency  towards  hysteresis  for  extreme 
conditions.  Unsteady  flow  effects  were  not  considered  in  the  current  study  since 
neither  of  the  typical  cases  developed  had  high  enough  flows  to  magnify  unsteady 
flow  effects  to  the  degree  demonstrated  by  the  Mississippi  River  study  (Westphal,  et 
al.,  1999). 
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3.3.3  DAMAGE-STAGE  UNCERTAINTY 

Sensitivity  to  damage-stage  uncertainty  was  examined  in  a  similar  manner  to  stage- 
flow  uncertainty.  A  normal  distribution  was  assumed  to  describe  the  uncertainty,  with 
the  standard  deviation  of  error  about  the  mean  variable.  The  standard  deviation  of 
error  was  assumed  to  increase  linearly  with  stage,  with  zero  uncertainty  at  (0,0).  Due 
to  the  intense  data  requirements  of  doing  a  rigorous  derivation  of  the  damage-stage 
function,  the  effect  of  large  variances  in  damages  were  examined  as  recommended  by 
Lindquist  (1995).  EAD  and  EAD  reduction  were  calculated  for  standard  deviation  of 
error  equal  to  5%,  10%  and  20%  of  the  mean  damage  value. 

Venkatesh  and  Hobbs  (1999)  assumed  a  uniform  uncertainty  distribution  for  a  study  to 
compare  the  effects  of  climate  uncertainty  with  damage  uncertainty  at  Lake  Erie.  The 
uniform  distribution  had  1/3  probability  that  the  damage-stage  curve  will  be  50% 
lower  than  its  expected  value,  1/3  probability  that  the  damage-stage  curve  will  be  50% 
higher  than  its  expected  value,  and  1/3  probability  that  the  damage-stage  curve  will  be 
equal  to  its  expected  value.  This  assumption  is  different  than  for  the  present  study  in 
that  the  present  study  assumes  that  enough  knowledge  is  present  that  the  best  estimate 
is  known  to  be  more  likely  than  surrounding  values. 

3.4  Uncertainty  Analysis 

An  uncertainty  analysis  was  performed  to  compare  the  relative  contributions  of 
uncertainty  in  each  of  the  primary  input  functions  to  uncertainty  in  the  resulting  EAD 
without  a  levee  and  with  50,  100  and  250-year  (0.02,  0.01  and  0.004  exceedance 
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probability)  levees.  Uncertainty  was  varied  in  each  input  function  according  to  the 
ranges  defined  in  the  Sensitivity  to  Uncertainty  section.  Uncertainty  in  the  output  was 
defined  as  the  difference  between  the  0.25  exceedance  probability  and  the  0.75 
exceedance  probability  EAD  (Goldman,  pers.  comm.,  2000). 

3.5  EAD  Reduction 

To  examine  the  sensitivities  of  EAD  reduction,  with-project  conditions  were  created. 
EAD  reduction  was  examined  for  levees  built  to  the  50-year,  100-year,  and  250-year 
(0.02,  0.01  and  0.004  exceedance  probability)  levels.  The  same  set  of  runs  was  done 
for  the  with-project  conditions  as  for  the  without-project  condition.  The  without- 
project  EAD  values  were  subtracted  from  the  with-project  EAD  values  to  get  EAD 
reduction,  an  estimate  of  benefits  expected  from  the  levee  project. 

3.6  Sensitivity  Index 

To  analyze  the  sensitivity  of  EAD,  EAD  reduction,  and  their  associated  uncertainties 
to  the  above-described  parameters  on  a  consistent  basis,  a  sensitivity  index  must  be 
calculated.  Simple  sensitivity  can  be  used  as  a  first  measure  of  uncertainty 
importance.  Mathematically,  simple  sensitivity  is  defined  as: 

Us(x,y)  =  [8y/5x]X0 
(Morgan  and  Henrion,  1 990) 

In  terms  of  flood  damage  analysis,  simple  sensitivity  can  be  written  as: 

Sxi(xi,x2,x3,y)  =  [8y/8xi]X]0,X2o,x3o 
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where  S  =  Simple  Sensitivity  of  EAD  to  xj 

X]  =  Hydrologic  Function;  X2  =  Hydraulic  Function;  X3  =  Economic  Function;  y  = 
EAD; 

xio  =  Hydrologic  Function  at  Base  Case 
X20  =  Hydraulic  Function  at  Base  Case 
X30  =  Economic  Function  at  Base  Case. 

This  definition  works  well  for  parameters  on  the  same  scale.  Simple  sensitivity  can  be 
used  to  examine  changes  in  input  function  sensitivity  over  a  range  of  values. 
However,  since  the  three  primary  relationships  in  flood  damage  analysis  are  on 
different  scales,  a  normalized  sensitivity  must  be  used  to  compare  the  sensitivities 
between  functions.  The  normalized  sensitivity  is  also  known  as  elasticity,  which  is 
essentially  the  percent  change  in  output  with  the  percent  change  in  the  uncertain 
parameter  or  input,  defined  as: 

UE(x,y)  =  [Sy/8x]X0  (x0/y0) 

(Morgan  and  Henrion,  1990) 

In  terms  of  flood  damage  analysis,  this  can  be  rewritten  as: 

Exi(X],x2,x3,y)  =  [8y/8xi]xio,x2o,x3o(xi0/y0) 
where  E  =  Elasticity  of  EAD  to  xj 

X]  =  Hydrologic  Function;  X2  =  Hydraulic  Function;  x3  =  Economic  Function;  y  = 
EAD; 

X]0  =  Hydrologic  Function  at  Base  Case 
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X20  =  Hydraulic  Function  at  Base  Case 

X30  =  Economic  Function  at  Base  Case;  y0  =  EAD  at  Base  Case. 

Simple  sensitivity  can  be  plotted  for  comparison  of  function  parameter  sensitivity  over 
the  range  of  values  examined.  Elasticity  can  be  plotted  for  comparison  of  EAD 
sensitivity  to  each  function.  The  same  plots  can  be  created  for  EAD  reduction  and  the 
variance  of  EAD. 

3.7  Case  Studies 

To  test  the  results  of  the  numerical  experiments  on  the  synthetic  data  sets,  sensitivity 
to  uncertainty  was  examined  for  two  sample  data  sets.  A  data  set  from  a  typical  large 
drainage  basin  on  the  Blue  River,  Missouri  and  a  data  set  from  a  typical  small 
drainage  basin  on  the  Chippewa  River,  Georgia  were  used  for  comparison.  Both  data 
sets  had  flow-exceedance  probability  curves  that  were  well  defined  by  an  LPIII 
distribution  and  previously  defined  stage-flow  and  damage-stage  curves. 
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Chapter  4 
Results 

The  following  graphs  summarize  the  results  of  the  Sensitivity  and  Uncertainty 
Analyses,  as  described  in  the  Approach  section.  Elasticity  (percent  change  in  EAD  or 
EAD  Reduction  divided  by  percent  change  in  parameter  value)  was  plotted  against 
parameter  value  so  the  sensitivities  could  be  compared  on  a  unitless  basis.  All 
numerical  results  are  included  in  Appendix  A.  A  vertical  line  was  drawn  on  each 
graph  to  delineate  the  base  case  value  against  which  elasticity  was  measured.  The 
sensitivity  analysis  results  are  presented  according  to  function  in  the  following  order: 

1 .  Flow-Exceedance  Probability  Mean  -  Figures  8-11 

2.  Flow-Exceedance  Probability  Standard  Deviation  -  Figures  12-14 

3.  Flow-Exceedance  Probability  Skew  -  Figures  15-17 

4.  Stage-Flow  Function  Average  Slope,  “B”  -  Figures  18-21 

5.  Stage-Flow  Function  Y-Intercept,  “C”  -  Figures  22  and  23 

6.  Damage-Stage  Function  Lower  Bound  of  Damage-Producing  Stages  -  Figures  24 
and  25 

7.  Damage-Stage  Function  Upper  Bound  of  Damage-Producing  Stages  -  Figures  26- 
28 

8.  Damage-Stage  Function  Inflection  Points,  A  and  B  -  Figure  29  -  32 
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4.1  Flow-Exceedance  Probability  Sensitivity  Analysis 

The  sensitivity  of  EAD  and  EAD  reduction  to  the  three  moments  of  the  flow- 
exceedance  probability  function  was  examined  by  computing  elasticity  relative  to  the 
large  and  small  base  cases.  The  elasticity  of  EAD  and  EAD  reduction  to  mean  flow 
for  the  large  basin  was  the  highest  for  the  first  few  thousand  cubic  feet  per  second 
below  the  base  case,  and  declined  as  the  flow  decreased.  The  small  basin  showed  a 
steadier  decline  in  elasticity  as  flow  decreased.  The  elasticity  of  EAD  and  EAD 
reduction  to  mean  flow  was  higher  for  the  small  basin  than  for  the  large  basin  over 
nearly  the  entire  flow  range  examined.  For  both  basins,  EAD  reduction  with  a  50-year 
(0.02  exceedance  probability)  levee  was  more  sensitive  to  mean  flow  than  EAD 
reduction  with  the  larger  levees.  As  the  levee  size  increases,  the  elasticity  of  EAD 
reduction  approaches  that  of  EAD  without  a  levee.  As  levee  size  increases,  the 
damage  prevented  approaches  the  total  damage  expected  without  a  project. 

While  the  elasticity  of  EAD  and  EAD  reduction  with  a  50-year  (0.02  exceedance 
probability)  levee  increased  steadily  with  mean  flow,  the  elasticity  of  EAD  reduction 
with  100  and  250-year  (0.01  and  0.004  exceedance  probability)  levees  decreased 
temporarily  at  a  mean  flow  of  around  24,000  cfs.  This  occurs  because  the  average 
steepness  of  the  damage-producing  part  of  the  stage-damage  curve  decreases  when 
mean  flow  reaches  24,000  cfs  with  a  100  or  250-year  (0.01  or  0.004  exceedance 
probability)  levee.  As  the  mean  flow  is  increased  to  24,500  cfs,  the  average  steepness 
of  the  damage-producing  part  of  the  stage-damage  curve  does  not  change  compared  to 
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what  it  was  at  a  mean  of  24,000  cfs,  so  the  elasticity  increases  again  due  to  the 
increase  in  flow.  The  results  are  summarized  in  Figures  8  and  9. 
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Elasticity  of  EAD  to  Mean  Flow  -  Small  Basin 
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Figure  9.  Elasticity  of  EAD  and  EAD  Reduction  with  Mean  Flow,  Small 
Basin 
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The  change  in  mean  flow  is  typically  propagated  through  the  EAD  calculation  through 
a  series  of  steps.  If  new  data  were  obtained  that  decreased  the  best  estimate  of  the 
mean  flow,  the  entire  flow-exceedance  probability  curve  would  be  lowered,  thus 
changing  the  flow  range  required  for  the  stage-flow  curve  to  be  integrated  with  the 
flow-exceedance  probability  curve.  To  isolate  the  mean  for  sensitivity  analysis,  only 
the  range  of  the  stage-flow  curve  was  changed  to  match  the  range  of  flows  for  each 
mean  flow.  However,  examination  of  data  sets  from  across  the  United  States  showed 
that  mean  flow  can  be  correlated  with  the  average  slope  of  the  rating  curve,  as  shown 
in  Figures  5a  and  5b.  Changing  the  average  slope  of  the  rating  curve  according  to  this 
correlation  with  mean  flow  changes  the  behavior  of  EAD,  EAD  reduction  and 
elasticity. 

The  mean  flow  -  stage-flow  curve  slope  correlation  shows  that  streams  with  lower 
mean  flows  tend  to  have  steeper  rating  curves.  A  steeper  rating  curve  tends  to 
increase  EAD,  while  a  lower  mean  flow  will  tend  to  decrease  EAD.  A  set  of  model 
runs  was  completed  to  demonstrate  the  interaction  of  these  tendencies.  The  resulting 
elasticity  changes  can  be  seen  in  Figures  10  and  11.  When  the  rating  curve  becomes 
steep  enough,  the  effect  of  steepening  the  curve  dominates  the  EAD  elasticity.  When 
the  steepness  of  the  rating  curve  dominates  EAD  elasticity,  the  elasticity  can  decrease. 
This  is  most  evident  in  the  small  basin  (Figures  9  and  11).  When  the  rating  curve  is 
flat  enough  the  mean  flow  dominates  the  EAD  elasticity,  increasing  its  values  at  the 
higher  flows.  This  is  most  evident  in  the  large  basin  (Figures  8  and  10). 
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Figure  11.  Elasticity  of  EAD  and  EAD  Reduction  to  Mean  Flow  with 
Adjusted  Rating  Curve,  Small  Basin 
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The  elasticity  of  EAD  and  EAD  reduction  with  standard  deviation  of  flow  was  the 
highest  for  the  first  few  cubic  feet  per  second  above  the  base  case,  and  decreased 
steadily  for  both  the  large  and  small  basins.  EAD  and  EAD  reduction  for  all  levee 
sizes  were  more  elastic  for  changes  in  standard  deviation  of  the  flow-exceedance 
probability  function  than  for  the  mean  of  the  flow  exceedance  probability  function 
over  the  ranges  examined.  As  with  the  mean  flow,  EAD  and  EAD  reduction  were 
more  elastic  for  the  large  basin  than  for  the  small  basin.  Also,  the  trend  of  decreasing 
elasticity  with  increasing  levee  size  was  opposite  of  the  previous  results.  These  results 
are  summarized  in  Figures  12  and  13. 
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Elasticity  of  EAD  to  Flow  Frequency  Std  Deviation  -  Large  Basin 
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Figure  12.  Elasticity  of  EAD  and  EAD  Reduction  with  Standard  Deviation 
of  Flow,  Large  Basin 


Elasticity  of  EAD  to  Flow  Frequency  Standard  Deviation  -  Small  Basin 
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Figure  13.  Elasticity  of  EAD  and  EAD  Reduction  with  Standard  Deviation 
of  Flow,  Small  Basin 


5 


Changes  in  the  standard  deviation  of  the  flow-exceedance  probability  function  are 
propagated  differently  than  changes  in  the  mean  flow.  As  the  standard  deviation  is 
increased,  the  slope  of  the  flow-exceedance  probability  function  increases.  This 
creates  two  opposing  trends,  as  illustrated  in  Figure  14. 

1 .  Less  frequent  (high  flow)  events  are  larger. 

2.  More  frequent  events  are  smaller. 

Construction  of  a  levee  eliminates  damages  from  the  more  frequent  decreased  lower 
flows,  but  the  damages  from  the  increased  higher  flows  remain.  The  net  effect  of 
increasing  the  standard  deviation  of  flow  is  to  increase  EAD  and  EAD  reduction,  but 
decrease  elasticity. 


Figure  14.  The  Effects  of  Increasing  Standard  Deviation  of  the  Flow* 
Exceedance  Probability  Function 
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The  elasticity  of  EAD  and  EAD  reduction  with  the  flow-exceedance  probability  skew 
follows  the  same  behavior  for  the  large  basin  as  for  the  small  basin.  Positive  skew 
values  increase  the  magnitude  of  the  highest  probability  and  lowest  probability  flows, 
while  decreasing  the  magnitude  of  the  middle  probability  flows  to  a  lesser  degree. 
Negative  skew  values  decrease  both  of  these  magnitudes,  while  increasing  the 
magnitude  of  the  middle  probability  flows  to  a  lesser  degree.  These  trends  can  be  seen 
in  Figure  15. 


The  Effects  of  Skew  on  the  Flow-Exceedance  Probability  Function 


Figure  15.  The  Effects  of  Skew  on  the  Flow-Exceedance  Probability 
Function 


Compared  to  the  mean  and  standard  deviation  of  the  flow-exceedance  probability 
function,  elasticity  decreases  very  slowly  with  increasing  skew  for  the  small  basin,  and 
even  more  slowly  for  the  large  basin.  For  skew,  elasticity  was  measured  against  a 
base  value  of  0.1  because  elasticity  cannot  be  computed  against  a  value  of  zero.  There 
are  ranges  of  negative  skew  over  which  the  elasticities  of  EAD  and  EAD  reduction 
remain  approximately  constant  or  increase  very  slightly.  This  occurs  when  increasing 
the  higher  flows  can  no  longer  affect  EAD  because  there  is  already  100%  damage.  In 
this  case,  EAD  is  affected  slightly  by  the  change  in  the  middle  probability  flows  in  the 
opposite  direction.  When  the  skew  coefficient  was  small  enough  to  shift  the  flow- 
exceedance  probability  curve  below  the  100-year  (0.01  exceedance  probability)  and 
250-year  (0.004  exceedance  probability)  levee  heights  that  were  sized  at  zero  skew  the 
EAD  became  approximately  zero.  Since  this  is  unrealistic,  the  flow-exceedance 
probability  curve  is  not  extended  beyond  the  0.001  exceedance  probability.  Over  all 
skew  ranges  examined,  both  the  large  and  small  basins  were  significantly  less 
sensitive  to  skew  than  to  mean  or  standard  deviation  of  the  flow-exceedance 
probability  function.  The  results  are  summarized  in  Figures  16  and  17. 


54 


Figure  17.  Elasticity  of  EAD  and  EAD  Reduction  with  Skew  of  Flow 
Exceedance  Probability  Function,  Small  Basin 
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4.2  Stage-Flow  Sensitivity  Analysis 

The  average  slope  of  the  stage-flow  function  is  defined  by  the  “b”  parameter  in  the 
general  quadratic  equation  used  to  define  the  stage- flow  curve,  as  discussed  in  Chapter 
3: 

S  =  aQ2  +  bQ  +  c  where  S  =  Stage  (ft)  and  Q  =  Flow  (cfs). 

The  elasticity  of  EAD  and  EAD  reduction  with  the  average  slope  of  the  stage-flow 
function  is  greatest  for  steeper  functions.  It  follows  that  the  elasticity  is  greatest  for 
smaller  basins.  In  general,  as  the  average  slope  increases,  EAD  increases  quickly.  As 
the  slope  of  the  stage-flow  function  increases,  higher  stages  become  associated  with 
the  same  flows,  and  thus  higher  damages  until  damages  reach  a  maximum  limit.  The 
elasticity  of  EAD  increases  initially  as  slope  increases  above  the  base  case,  then 
decreases  as  EAD  approaches  total  damage.  The  elasticity  of  EAD  reduction 
decreases  steadily  with  increasing  slope  above  the  base  case.  As  the  average  slope 
becomes  flatter  than  the  base  case,  elasticity  of  EAD  reduction  quickly  approaches 
that  of  EAD.  The  results  are  summarized  in  Figures  18  and  19.  The  range  in  stage- 
flow  function  slopes  examined  can  be  seen  in  Figures  20  and  21 . 
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Figure  18.  Elasticity  of  EAD  and  EAD  Reduction  with  Stage-Flow 
Function  Average  Slope,  Large  Basin 


Figure  19.  Elasticity  of  EAD  and  EAD  Reduction  with  Stage-Flow 
Function  Average  Slope,  Small  Basin 
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Stage  Flow  Function  Average  Sample  Slope  Variation,  Large 


Figure  20.  Variation  in  Stage-Flow  Function  Slope,  Large  Basin 


Stage  Flow  Function  Average  Sample  Slope  Variation,  Small 


Figure  21.  Variation  in  Stage-Flow  Function  Slope,  Small  Basin 
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The  “C”  parameter  in  the  synthetic  stage-flow  functions  represents  the  y-intercept  of 
the  curve  as  given  in  the  general  quadratic  equation: 

S  =  aQ2  +  bQ  +  c  where  S  =  Stage  (ft)  and  Q  =  Flow  (cfs). 

Increasing  “C”  increases  EAD  and  EAD  reduction.  For  the  large  basin,  elasticity 
increased  with  increasing  “C”.  For  the  small  basin,  elasticity  increased  with 
increasing  “C”,  but  began  to  decrease  as  the  maximum  damage  was  approached.  In 
general,  the  elasticity  of  EAD  and  EAD  reduction  to  the  “C”  parameter  was  similar  for 
each  basin  size.  The  results  are  shown  in  Figures  22  and  23.  The  change  in  “C”  value 
is  propagated  the  same  way  as  the  change  in  function  slope.  Increasing  “C”  increases 
the  stage  for  each  flow,  thereby  increasing  the  damages  until  a  maximum  limit  is 
reached. 
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Elasticity  of  EAD  to  Stage  Flow  Function  C  Parameter  -  Large  Basin 


Figure  22.  Elasticity  of  EAD  and  EAD  Reduction  with  Stage-Flow 
Function  “C”  Parameter,  Large  Basin 


Elasticity  of  EAD  to  Stage  Flow  Function  C  Parameter  -  Small  Basin 


Figure  23.  Elasticity  of  EAD  and  EAD  Reduction  with  Stage-Flow 
Function  “C”  Parameter,  Small  Basin 


4.3  Damage-Stage  Sensitivity  Analysis 

Elasticity  was  calculated  for  a  range  of  lower  bounds  of  the  damage-producing  stages. 
The  lower  bound  is  the  stage  at  which  damage  begins.  As  the  lower  bound  of  damage- 
producing  stages  was  increased,  EAD  and  EAD  reduction  decreased.  Elasticity  also 
decreased  steadily  with  an  increasing  lower  bound  for  both  the  large  and  small  basins. 
As  the  lower  bound  increased,  the  levees  prevented  a  greater  proportion  of  damages 
from  occurring.  The  elasticity  of  EAD  and  EAD  reduction  was  more  than  twice  as 
high  for  the  small  basin  as  for  the  large  basin  for  the  values  examined.  Increasing  the 
lower  bound  effectively  eliminates  damage  from  the  most  common  events,  while 
changing  the  less  common  events  relatively  little.  The  results  are  summarized  in 
Figures  24  and  25. 
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Figure  24.  Elasticity  of  EAD  and  EAD  Reduction  with  Damage-Stage 
Function  Lower  Bound  of  Damages,  Large  Basin 


Figure  25.  Elasticity  of  EAD  and  EAD  Reduction  with  Damage-Stage 
Function  Lower  Bound  of  Damages,  Small  Basin 
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Elasticity  was  calculated  for  each  basin  with  different  upper  bound  stages.  The  upper 
bound  stage  is  the  stage  above  which  maximum  damage  occurs.  As  the  upper  bound 
was  increased,  EAD  and  EAD  reduction  decreased.  Elasticity  also  decreased  steadily 
for  the  large  basin,  but  at  different  rates  depending  on  how  high  the  levee  was  for  the 
small  basin  (See  Figures  27  and  28).  Considering  that  damage  is  assumed  to  only 
occur  above  the  top  of  levee  stage,  the  average  relative  steepness  of  the  damage- 
producing  part  of  the  curve  varies  with  different  levee  stages  (See  Figure  26).  For 
instance,  without  a  levee  the  curve  for  the  upper  bound  equal  to  6.2  is  steeper  than  the 
curve  with  7.2.  However,  when  a  50-year  (0.02  exceedance  probability)  levee  is 
added,  the  damage-producing  part  of  the  curve  for  the  upper  bound  equal  to  6.2 
becomes  flatter  than  the  curve  with  7.2.  When  a  100-year  (0.01  exceedance 
probability)  levee  is  added,  the  curve  with  6.2  remains  flatter,  but  when  a  250-year 
(0.004  exceedance  probability)  levee  is  added  it  becomes  steeper  again.  These 
changes  are  mirrored  in  the  elasticity  curve  for  the  small  basin  (Figure  28). 
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Stage  (ft) 


Figure  26.  Variation  of  Damage-Stage  Curve  Upper  Bound 
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Figure  27.  Elasticity  of  EAD  and  EAD  Reduction  with  Damage-Stage 
Function  Upper  Bound  of  Damage,  Large  Basin 
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Figure  28.  Elasticity  of  EAD  and  EAD  Reduction  with  Damage-Stage 
Function  Upper  Bound  of  Damage,  Small  Basin 
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Elasticity  was  calculated  for  variation  in  each  inflection  point  of  the  damage-stage 
function  independently.  Since  only  three  values  were  examined  for  the  first  (lower) 
and  second  (higher)  inflection  points  apiece,  the  elasticity  values  were  graphed  as  a 
maximum  and  minimum  value.  In  general,  as  the  first  inflection  point  “A”  increased, 
EAD  decreased  for  both  the  large  and  small  basins.  As  the  second  inflection  point 
“B”  increased,  EAD  decreased  for  both  the  large  and  small  basins.  EAD  had  a  higher 
maximum  elasticity  for  the  small  basin  than  for  the  large  basin  for  both  inflection 
points,  but  a  lower  minimum  elasticity  for  both  inflection  points. 

EAD  and  EAD  reduction  had  a  higher  elasticity  in  both  basins  to  the  first  inflection 
point  than  to  the  second.  The  propagation  of  changing  an  inflection  point  in  the  stage 
damage  function  is  similar  to  changing  the  lower  bound  of  damage-producing  stage. 
In  general,  the  larger  levee  scenarios  had  a  lower  maximum  elasticity  of  EAD 
reduction,  but  a  higher  minimum  elasticity.  This  can  be  seen  in  Figures  29  —  32. 
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Figure  29.  Elasticity  of  EAD  with  Damage-Stage  Function  Inflection  Point 
Locations  for  Large  and  Small  Basins 
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Figure  30.  Elasticity  of  50-yr  Levee  EAD  Reduction  with  Damage-Stage 
Function  Inflection  Point  Locations  for  Large  and  Small  Basins 
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Figure  31.  Elasticity  of  100-yr  Levee  EAD  Reduction  with  Damage-Stage 
Function  Inflection  Point  Locations  for  Large  and  Small  Basins 
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Figure  32.  Elasticity  of  250-yr  Levee  EAD  Reduction  with  Damage-Stage 
Function  Inflection  Point  Locations  for  Large  and  Small  Basins 
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4.4  Sensitivity  to  Uncertainty 

Based  on  the  results  of  HEC-FDA  for  a  range  of  equivalent  record  lengths,  elasticity 
was  calculated  for  each  record  length  to  compare  results  on  a  unitless  basis.  Each 
individual  model  run  output  is  included  in  Appendix  A.  Elasticity  values  are 
summarized  in  Figure  33  for  the  small/steep  and  large/flat  basins.  Likewise,  elasticity 
was  calculated  for  each  standard  deviation  of  error  about  the  flow-stage  curve,  with 
each  individual  model  run  output  included  in  Appendix  A.  Elasticity  values  for  flow- 
stage  uncertainty  are  summarized  in  Figure  34.  Since  the  mean  EAD  is  not  sensitive 
to  stage-damage  uncertainty  (Lindquist,  1995),  no  plot  was  produced  for  its  elasticity. 
However,  model  runs  were  still  completed  to  verify  this  result  and  are  included  in 
Appendix  A. 
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Figure  33.  EAD  Elasticity  to  Equivalent  Record  Length 


Figure  34.  EAD  Elasticity  to  Stage-Flow  Error 
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Over  the  range  of  values  examined,  EAD  had  the  highest  elasticity  with  flow  record 
length.  This  was  always  the  case  for  the  large  basin,  and  usually  the  case  for  the  small 
basin.  EAD  had  the  highest  elasticity  with  stage-flow  standard  deviation  of  error  for 
the  small  basin  when  the  record  length  was  long  and  the  standard  deviation  of  stage- 
flow  error  was  large. 


The  maximum  and  minimum  elasticity  of  EAD  to  each  function  uncertainty  was 
calculated  from  the  range  of  results.  These  values  can  be  compared  for  the  large  and 
small  basins  in  Figures  35. 


Figure  35.  Uncertainty  Elasticity  Comparison  —EAD 
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Based  on  the  results  from  varying  record  length  with  50,  100  and  250-year  (0.02,  0.01 
and  0.004  exceedance  probability)  levees,  elasticity  to  equivalent  record  length  was 
calculated  for  EAD  reduction,  the  difference  between  EAD  without  and  with  each 
levee.  Elasticity  of  EAD  reduction  to  equivalent  record  length  was  plotted  for  each 
levee  in  each  synthetic  basin  in  Figure  36.  Elasticity  was  also  computed  for  the 
standard  deviation  of  error  of  the  stage-flow  function,  and  plotted  in  Figure  37.  The 
elasticity  of  EAD  reduction  generally  decreases  with  increasing  record  length,  but 
increases  slightly  in  some  cases  depending  on  the  slope  of  the  damage-stage  curve. 
The  changing  of  the  average  slope  of  the  damage-stage  curve  can  be  seen  in  Figure  26. 
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Elasticity  of  EAD  Reduction  Versus  Record  Length  for 
50,  100,  and  250-yr  Levees 


Figure  36.  EAD  Reduction  Elasticity  to  Equivalent  Record  Length 


The  behavior  of  the  elasticity  of  EAD  reduction  was  similar  to  that  of  EAD,  but  on  a 
smaller  scale.  EAD  reduction  had  the  largest  elasticity  with  respect  to  equivalent 
record  length  over  the  entire  range  for  the  large  basin,  and  over  most  of  the  range  for 
the  small  basin.  However,  for  long  record  lengths  and  large  standard  deviations  of 
error  for  the  flow  stage  function,  EAD  reduction  was  more  elastic  with  respect  to 
stage-flow  function  uncertainty.  The  ranges  of  elasticity  of  EAD  reduction  for  the  50, 
100  and  250-year  (0.02,  0.01  and  0.004  exceedance  probability)  levees  are 
summarized  in  Figures  38  through  40. 
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Figure  38.  Uncertainty  Elasticity  Comparison  -EAD  Reduction  with  50- 
Year  (0.02  Exceedance  Probability)  Levee 
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Figure  39.  Uncertainty  Elasticity  Comparison  -EAD  Reduction  with  100 
Year  (0.01  Exceedance  Probability)  Levee 
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Figure  40.  Uncertainty  Elasticity  Comparison  -EAD  Reduction  with  250 
Year  (0.004  Exceedance  Probability)  Levee 
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4.5  Uncertainty  Analysis 

An  uncertainty  analysis  was  completed  to  examine  the  contribution  of  uncertainty  in 
each  primary  function  to  uncertainty  in  EAD  for  no  levee  and  with  50,  100  and  250- 
year  (0.02,  0.01  and  0.004  exceedance  probability)  levees.  To  provide  a  basis  for 
comparison,  uncertainty  in  EAD  was  defined  as  the  difference  between  the  0.25 
exceedance  and  0.75  exceedance  quantiles  in  the  output  distribution.  The  numerical 
results  are  included  in  Appendix  A  and  graphed  in  Figures  41, 42,  and  43  for  the  flow- 
exceedance  probability,  stage-flow  and  damage-stage  functions,  respectively. 


Difference  Between  25%  Exceedance  and  75%  Exceedance  Values  of  EAD  With  Record 

Length 


Figure  41.  Flow-Exceedance  Probability  Uncertainty  Analysis 
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Difference  Between  25%  Exceedance  and  75%  Exceedance  Values  of  EAD  With  Standard 
Deviation  of  Stage  Flow  Error 


Figure  42.  Stage-Flow  Uncertainty  Analysis 
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The  relative  contribution  of  uncertainty  in  each  primary  function  to  uncertainty  in 
EAD  depends  upon  in  which  ranges  the  individual  uncertainties  lie.  The  relative 
contribution  of  uncertainty  increases  as  record  length  decreases  or  the  standard 
deviation  of  error  of  the  stage-flow  function  or  damage-stage  function  increases.  In 
general,  the  uncertainty  contribution  for  the  small  basin  is  larger  than  the  uncertainty 
contribution  for  the  large  basin  with  the  same  input  uncertainty.  This  difference  is 
smaller  for  the  flow-exceedance  probability  and  damage-stage  functions  and  more 
pronounced  for  the  stage-flow  function.  The  contribution  of  each  function  to 
uncertainty  in  EAD  decreases  with  increasing  levee  size.  The  uncertainty  contribution 
from  the  flow-exceedance  probability  function  begins  to  level  off  for  record  lengths 
greater  than  about  60  years  without  a  levee,  and  greater  than  about  50  years  with  a 
levee.  The  uncertainty  contribution  from  the  stage-flow  function  increases  more 
steadily  with  standard  deviation  of  error,  as  does  the  contribution  from  the  damage- 
stage  function. 

Over  most  of  the  ranges  examined,  stage-flow  uncertainty  appears  to  have  the  greatest 
contribution  to  uncertainty  in  EAD  for  the  small  basin,  whereas  flow-exceedance 
probability  uncertainty  appears  to  have  the  greatest  contribution  to  uncertainty  in  EAD 
for  the  large  basin.  While  the  relative  strength  of  each  contribution  can  be  inferred  by 
isolating  uncertainties  in  this  manner,  the  interaction  between  uncertainties  is  not  as 
easily  defined. 
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4.6  Case  Study  Results 

The  results  of  the  Chippewa  River  and  Blue  River  case  studies  verified  the  results  of 
the  sensitivity  to  uncertainty  experiments  conducted  on  the  synthetic  basins.  Slight 
differences  in  elasticity  can  be  attributed  mainly  to  the  different  shape  of  the  damage- 
stage  curves  in  the  case  studies.  The  synthetic  basin  experiments  showed  that  EAD  is 
sensitive  to  the  location  of  the  upper  and  lower  bounds  of  the  damage-stage  functions. 
In  general,  the  case  studies  verified  that  the  large  basin  is  more  sensitive  to  flow- 
exceedance  probability  uncertainty  due  to  record  length  than  to  stage-flow  uncertainty. 
The  small  basin  is  usually  more  sensitive  to  flow-exceedance  probability  uncertainty, 
except  for  situations  of  long  record  length  and  high  stage-flow  uncertainty.  The  large 
basin  is  always  more  sensitive  to  record  length  than  the  small  basin,  and  the  small 
basin  is  always  more  sensitive  to  stage-flow  uncertainty  than  the  large  basin  for  the 
ranges  of  values  examined.  The  case  study  results  can  be  seen  in  Figures  44  and  45, 
and  compared  to  the  results  of  the  synthetic  basins  in  Figures  33  and  34. 


Elasticity  of  EAD  Versus  Record  Length  -  Blue  River  and  Chippewa  River  Case  Studies 


Figure  44.  Elasticity  of  EAD  Versus  Record  Length  -  Case  Studies 
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Elasticity  of  EAD  Versus  Stage-Flow  Uncertainty 
Blue  River  and  Chippewa  River  Case  Studies 
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Figure  45.  Elasticity  of  EAD  Versus  Stage-Flow  Uncertainty  -  Case 
Studies 
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The  sensitivity  of  EAD  reduction  to  uncertainty  for  the  case  studies  followed  generally 
the  same  trends  as  for  the  synthetic  basins,  with  differences  due  to  the  shape  of  the 
damage-stage  function.  The  trends  of  EAD  elasticity  differed  some  from  the  trends  of 
EAD,  depending  on  the  height  of  the  levee.  While  both  synthetic  damage-stage 
functions  were  piecewise  linear  approximations  to  an  S-shaped  curve  that  started  and 
ended  at  the  same  frequencies  as  one  another,  the  case  studies  had  different  shapes  and 
started  and  ended  at  different  frequencies.  The  case  study  damage-stage  curves  are 
shown  in  Figures  46  and  47. 


Chippewa  River  Damage-Stage  Curve 


Figure  46.  Chippewa  River  Case  Study  Damage-Stage  Curve 
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Blue  River  Damage-Stage  Function 


Figure  47.  Blue  River  Case  Study  Damage-Stage  Curve 

While  the  synthetic  damage-stage  functions  both  begin  at  the  mean  stage  associated 
with  the  2-year  (0.5  exceedance  probability)  event  and  maximize  at  the  500-year 
(0.002  exceedance  probability)  event,  the  case  study  functions  cover  different  ranges. 
The  Chippewa  River  damage-stage  function  begins  at  the  mean  stage  associated  with 
the  1.4-year  (0.7  exceedance  probability)  event  and  maximizes  at  the  mean  stage 
associated  with  the  500-year  event  (0.002  exceedance  probability).  The  Blue  River 
damage-stage  function  begins  at  the  mean  stage  associated  approximately  with  the  1  - 
year  event  (0.999  exceedance  probability)  and  maximizes  at  the  mean  stage  associated 
with  the  500-year  (0.002  exceedance  probability)  event.  These  different  beginning 
and  ending  points  of  the  damage-stage  curves  affect  EAD,  EAD  reduction  and  their 
elasticity  as  discussed  in  the  sensitivity  analysis  on  the  lower  and  upper  bounds  of  the 
synthetic  damage-stage  functions.  The  combined  effects  of  different  shapes  and 
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ranges  of  the  damage-stage  curves  are  the  main  causes  of  the  irregularities  and 
differences  in  the  elasticity  trends,  as  shown  in  Figures  48  through  51.  Case  study 
results  in  Figures  48  through  51  can  be  compared  with  the  synthetic  case  results  in 
Figures  36  and  37. 
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Figure  48.  Elasticity  of  EAD  Reduction  Versus  Record  Length,  Chippewa 
River  (Small  Basin) 


Figure  49.  Elasticity  of  EAD  Reduction  Versus  Stage-Flow  Uncertainty, 
Chippewa  River  (Small  Basin) 
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Elasticity  of  EAD  Versus  Record  Length  for  50,  100,  and  250-yr  Levees 
Blue  River  Case  Study 


- EAD  Reduction,  50yr  Levee 

—  A  —  EAD  Reduction,  lOOyr  Levee 
— *  -  EAD  Reduction,  250yr  Levee 
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Figure  50.  Elasticity  of  EAD  Reduction  Versus  Record  Length  -  Blue 
River  (Large  Basin) 


Elasticity  of  EAD  Reduction  Versus  Stage-Flow  Uncertainty  for  50, 100,  and  250-yr  Levees 
Blue  River  Case  Study 


Figure  51.  Elasticity  of  EAD  Reduction  Versus  Stage-Flow  Uncertainty - 
Blue  River  (Large  Basin) 
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Chapter  5 
Conclusions 


The  results  of  the  numerical  experiments  showed  the  relative  effect  of  each  variable  on 
EAD  and  EAD  reduction,  with  some  differences  between  the  large  basin  and  the  small 
basin.  In  summary,  EAD  and  EAD  reduction  were  sensitive  to  the  following  variables 
in  the  order  of  increasing  importance  given  in  Tables  7  and  8. 


Table  7.  Large  Basin  Order  of  Elasticity 


EAD 

EAD  Reduction 

1 .  Flow-Exceedance  Probability  Skew 

Flow-Exceedance  Probability  Skew 

2.  Damage-Stage  Inflection  Pt  “B” 

Damage-Stage  Inflection  Pt  “B” 

3.  Damage-Stage  Inflection  Pt  “A” 

Damage- Stage  Inflection  Pt  “A” 

4.  Flow-Exceedance  Probability  S.  Dev. 

Flow-Exceedance  Probability  Mean 

5.  Flow-Exceedance  Probability  Mean 

Flow-Exceedance  Probability  S.  Dev. 

6.  Stage-Flow  Function  “C”  Parameter 

Stage-Flow  Function  “C”  Parameter 

7.  Damage- Stage  Upper  Bound 

Damage- Stage  Upper  Bound 

8.  Damage- Stage  Lower  Bound 

Damage-Stage  Lower  Bound 

9.  Stage-Flow  Function  Slope,  “B” 

Stage-Flow  Function  Slope,  “B” 
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Table  8.  Small  Basin  Order  of  Elasticity 


EAD 

EAD  Reduction 

1 .  Flow-Exceedance  Probability  Skew 

Flow-Exceedance  Probability  Skew 

2.  Damage-Stage  Inflection  Pt  “B” 

Damage-Stage  Inflection  Pt  “B” 

3.  Flow-Exceedance  Probability  Mean 

Flow-Exceedance  Probability  Mean 

4.  Damage-Stage  Inflection  Pt  “A” 

Damage-Stage  Inflection  Pt  “A” 

5.  Flow-Exceedance  Probability  S.  Dev. 

Flow-Exceedance  Probability  S.  Dev. 

6.  Stage-Flow  Function  “C”  Parameter 

Stage-Flow  Function  “C”  Parameter 

7.  Damage-Stage  Upper  Bound 

Stage-Flow  Function  Slope,  “B” 

8.  Stage-Flow  Function  Slope,  “B” 

Damage-Stage  Upper  Bound 

9.  Damage-Stage  Lower  Bound 

Damage-Stage  Lower  Bound 

The  range  of  EAD  elasticity  calculated  for  the  sensitivity  analysis  can  be  seen  in 
Figure  52  and  for  sensitivity  to  uncertainty  in  Figure  53. 
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Figure  52.  Sensitivity  Analysis  Summary 
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As  for  the  first  set  of  numerical  experiments,  the  sensitivity  of  EAD  and  EAD 
reduction  to  each  type  of  uncertainty  depended  on  the  size  of  the  basin  and  initial 
value  of  uncertainty.  In  general,  EAD  and  EAD  reduction  are  most  sensitive  to  flow- 
exceedance  probability  uncertainty,  followed  by  stage-flow  uncertainty  and  damage- 
stage  uncertainty.  Elowever,  EAD  and  EAD  reduction  become  more  sensitive  to 
stage-flow  uncertainty  for  a  small  basin  with  a  record  length  of  less  than  about  55 
years,  and  a  standard  deviation  of  stage-flow  error  of  greater  than  about  1.5  feet  (See 
Figures  39,  40,  42  and  43).  These  results  are  consistent  with  the  uncertainty  analysis 
results,  which  show  that  stage-flow  uncertainty  has  the  greatest  contribution  to 
uncertainty  in  EAD  for  the  small  basin,  whereas  the  flow-exceedance  probability 
uncertainty  has  the  greatest  contribution  to  uncertainty  in  EAD  for  the  large  basin 
under  these  conditions  (See  Figures  47  through  49). 

In  general,  the  case  studies  verified  relative  sensitivities  to  uncertainty  for  sample  data 
sets.  However,  the  case  studies  also  showed  that  the  relative  sensitivities  change  for 
different  shaped  damage-stage  curves.  In  the  case  studies  used,  EAD  for  the  small 
basin  (Chippewa  River)  was  more  sensitive  to  flow-exceedance  probability 
uncertainty  over  a  larger  range  of  record  lengths  than  EAD  for  the  synthetic  small 
basin.  As  with  the  synthetic  basins,  EAD  for  the  large  basin  (Blue  River)  was  more 
sensitive  than  the  small  basin  (Chippewa  River)  to  flow-exceedance  probability 
uncertainty,  and  EAD  for  the  small  basin  was  more  sensitive  than  EAD  for  the  large 
basin  to  stage-flow  uncertainty. 
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Based  on  the  results  summarized  above,  several  trends  emerged  that  could  be  useful  in 
federal  investment  decision-making.  Variables  with  the  highest  elasticity  are  those 
that  have  the  most  potential  to  improve  estimates  of  EAD  and  EAD  reduction.  It  may 
be  worth  the  extra  expenditure  to  refine  the  best  estimate  of  a  particular  function  or  to 
reduce  the  uncertainty  if  the  cost  is  small  compared  to  the  benefits  received. 

5.1  Recommendations  for  Investment 

The  results  of  the  sensitivity  analyses  can  be  used  to  help  make  investment  decisions 
to  improve  estimates  of  EAD  and  EAD  reduction.  For  instance,  if  there  is  a  small  data 
set  (i.e.,  less  than  10-year  record)  and  funding  to  improve  the  estimate  of  EAD 
reduction,  it  may  be  wise  to  wait  another  year  or  more  to  reduce  the  hydrologic 
uncertainty.  The  expected  benefits  of  waiting  a  year  would  be  affected  by  a 
combination  of  the  reduction  in  hydrologic  uncertainty  and  a  slight  change  in  the  flow- 
exceedance  probability  function  parameters  (Davis,  et.  al.,  1972).  The  benefits  could 
be  estimated  as  the  difference  between  the  EAD  with  the  current  data  minus  the  EAD 
with  one  more  year  of  data.  The  expected  costs  of  waiting  another  year  are  a 
combination  of  any  damage  suffered  during  that  year  minus  the  annualized  cost 
avoided  by  deferring  construction  for  one  year  (Davis,  et.  ah,  1972). 

Along  the  same  lines,  it  might  be  wiser  to  refine  the  estimate  of  the  stage-flow  curve 
or  the  damage-stage  relationship  by  reducing  the  uncertainty  and  refining  the  most 
likely  curves.  The  benefits  of  improving  the  estimate  of  EAD  in  this  manner  can  be 
assessed  using  the  value  of  additional  information,  as  described  in  Chapter  2.  The 
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sensitivity  analyses  demonstrate  over  which  range  of  function  values  the  investor 
would  get  the  greatest  improvement  in  EAD  and  EAD  reduction  values,  and  hence  the 
most  benefit.  It  is  most  important  to  get  an  accurate  estimate  of  the  variables  with  the 
highest  elasticity,  such  as  stage-flow  function  shape,  and  least  important  to  get  an 
accurate  estimate  of  the  variables  with  the  lowest  elasticity,  such  as  flow-exceedance 
probability  skew. 

5.2  Regional  Differences 

Regional  differences  in  EAD  and  EAD  reduction  sensitivity  to  primary  input  function 
parameters  occur,  especially  if  the  LPIII  distribution  fits  regional  flow  data  poorly. 
Graphical  flow-exceedance  probability  curves  or  curves  better  defined  by  distributions 
other  than  the  LPIII  may  affect  the  elasticity  of  EAD  and  EAD  reduction  differently, 
and  affect  the  recommendations  for  investment.  Reservoir  regulation  and  flow 
augmentations  are  some  factors  that  can  contribute  to  these  differences. 

Regional  differences  in  EAD  and  EAD  reduction  sensitivity  to  primary  input  function 
parameters  also  occur  due  to  differences  in  the  shape  of  the  stage-flow  curve.  For 
instance,  water  in  a  narrow  channel  and  floodplain  will  tend  to  rise  faster  than  in  a 
wide  channel  and  floodplain  with  the  same  flow  capacity.  Whether  streamflow  is 
generated  by  short,  strong  bursts  of  rainfall  or  gradual  melting  of  snowpack  can  affect 
the  shape  of  the  channel  and  stage-flow  curve  significantly.  Likewise,  regional 
development  patterns  and  land  use  characteristics  are  the  main  contributors  to 
differences  in  stage-damage  curves  and  the  sensitivity  of  EAD  to  changes  in  the  curve. 
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5.3  Further  Study 

In  continuation  of  the  work  completed  for  this  study,  additional  numerical  experiments 
are  recommended  to  further  develop  Risk-based  Analysis.  Similar  tests  could  be 
completed  on  smaller  and  larger  basins,  different  shaped  (or  differently  defined) 
damage  curves,  or  different  flow-exceedance  probability  distributions.  Further  testing 
is  also  recommended  to  test  the  sensitivity  of  levee  reliability  measures  in  HEC-FDA. 
The  study  could  be  extended  to  include  uncertainty  with  flow  regulation. 
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Output  File  -  Sensitivity  Analysis 


CD 

Q_ 

O 

0 


CD 

E 


0) 


$ 

o 

M— 

I 

CD 
03 
CD 
-*— • 

CO 


"O 

c 

=3 

o 

-Q 

i_ 

0 

CL 

CL 

zs 


> 

0 

0 


CD 

E 

CD  » 

S>  0 

0  03 

—  0 

5  W 
o  i 
•*r  o) 
0  0 

CD  E 

S  § 
co  Q 


m  o 
co  _J 

LL  LL 

co  co 


z> 

CO 

CO 

Q 


DO 

CO 

CO 

o 


> 

0 

“O 

C  “O 
0  CO 

0  i 
^  c 
-  w 
c  0 
'co 

0  = 


0 


0 

_  E 

03  0 
</)  CT 

©  -1  2 
E  av 
#  2  5 
z  V  O 
c  £  u- 

3  O 

0£  ll 

75  " 

M  ^  co 
a  ~J  co 

>  LL.  LL 
H  LL  LL 


13  £ 

C  O 

o  o 

E  © 

^  m  c 
<D  g  - 

0  o  ^ 
.E  0"  0' 

0  CD  CD 

5  te  te 

a>  ®"  © 
fc:  o)  o) 
0  0  0 

-  0  0 

©  ©  © 
,s—  CD  CD 
0  0 

E  E 

0  0 


$ 
o 

ul  Q  Q 


S  =!  5 

U.  CO  CO 

LL  Q  Q 


CO 


0 

E 

0 

2: 

c 

ZD 

CZ 


O 

LO 

CO 

N- 

CO 

CM 

03 

CM 

CO 

CM 

CD 

CM 

00 

LO 

T— 

03 

CO 

LU 

00 

CO 

c\i 

CM 

c\i 

T- 

CM 

■ 

O 

O 

:: 

00 

CM 

00 

It* 

03 

CD 

CM 

LO 

CM 

LU 

LO 

CM 

LU 

M- 

''si- 

CO 

00 

CM 

CM 

CM 

O 

O 

0 

0 

0 

O 

O 

O 

LO 

O 

0 

0 

0 

O 

O 

O 

C 

X 

CM 

■ 

O 

0 

0 

0 

O 

O 

O 

5 

0 

CO 

■ 

d 

d 

d 

d 

d 

d 

d 

2 

O 

0 

■ 

0 

T— 

CO 

M- 

03 

CO 

T“ 

T” 

O 

■ 

0 

CD 

T— 

co 

T- 

03 

03 

O 

LU 

00 

00 

CM 

CM 

c\i 

T— 

T- 

O 

0 

03 

T— 

h- 

T — 

O 

0 

xt* 

00 

00 

00 

CM 

CM 

CM 

LU 

LU 

O 

O 

0 

0 

0 

O 

O 

O 

O 

O 

0 

0 

0 

O 

O 

O 

C 

A 

CD 

O 

0 

0 

0 

O 

O 

O 

CO 

■ 

d 

d 

d 

d 

d 

d 

d 

0 

■ 

CD 

00 

Tt 

LO 

CD 

0 

O 

LO 

h- 

T“ 

h- 

T— 

03 

03 

N- 

LU 

co 

00 

06 

CM 

CM 

CM 

T- 

00 

■ 

CD 

00 

0 

CD 

00 

CO 

O 

LO 

LU 

O 

LO 

CO 

00 

00 

CM 

CM 

CM 

T — 

UJ 

0 

0 

0 

O 

O 

O 

0 

O 

0 

0 

0 

O 

O 

O 

0 

C 

X 

0 

LO 

0 

0 

0 

O 

O 

O 

0 

CO 

d 

d 

d 

d 

d 

d 

d 

00 

03 

00 

T — 

00 

0 

CM 

CM 

00 

CM 

h- 

00 

03 

LO 

CM 

03 

03 

CM 

LU 

06 

00 

CM 

CM 

CM 

T- 

| 

CM 

CD 

CD 

Is- 

LU 

LU 

LO 

00 

00 

CM 

CM 

C 

X 

O 

0 

0 

0 

0 

O 

O 

0 

O 

0 

0 

0 

0 

O 

O 

O 

0 

0 

0 

0 

O 

O 

JO 

■ 

d 

d 

d 

d 

d 

d 

d 

O 

LO 

O 

0 

0 

0 

0 

0 

0 

0 

CM 

M" 

CM 

0 

0 

0 

0 

0 

0 

O 

O 

O 

0 

0 

0 

0 

0 

0 

< 

O 

O 

0 

0 

0 

0 

0 

0 

LU 

d 

d 

d 

d 

d 

d 

d 

d 

O 

O 

CD 

CM 

10 

0 

0 

0 

0 

0 

T- 

03 

LO 

CM 

0 

0 

0 

0 

0 

O 

O 

O 

0 

0 

0 

0 

0 

0 

< 

O 

O 

0 

0 

0 

0 

0 

0 

LU 

d 

d 

d 

d 

d 

d 

d 

d 

O 

03 

CD 

CD 

LO 

0 

0 

0 

0 

LO 

00 

0 

LO 

CM 

0 

0 

0 

0 

O 

0 

O 

0 

0 

0 

0 

< 

0 

0 

0 

O 

0 

0 

0 

0 

LU 

d 

d 

d 

d 

d 

d 

d 

d 

"k 

LO 

00 

LO 

LO 

00 

03 

00 

LO 

03 

03 

co 

00 

a 

00 

LO 

00 

CM 

T — 

0 

0 

0 

< 

0 

0 

0 

O 

O 

0 

0 

0 

LU 

d 

d 

d 

d 

d 

d 

d 

d 

&_ 

0 

0 

E 

CD 

l— 

0 

0 

D 

0 

0 

0 

0 

0 

0 

CL 

0 

0 

0 

LO 

0 

LO 

0 

LO 

<1 

> 

0 

LO 

CM 

;  cm 

00 

00 

L. 

0 

0 

E 

0 

CD 

l_ 

JD 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

LO 

0 

LO 

0 

LO 

0 

LO 

CL 

> 

N- 

CD 

CD 

LO 

10 

00 

100 


Output  File  -Sensitivity  Analysis 


CD 

CL 

O 


c 

“O 

c 

Cl 

c 

O 

o 

o 

u 

0 

q= 

o 

-Q 

c 

L_ 

C 

33 

0 

■*“ 

CL 

“O 

M— 

Cl 

c 

o 

O 

CNI 

> 

_ r 

_ r 

_0 

"ro 

*0 

CD 

E 

E 

CD 

0 

0 

CD 

i_ 

0~ 

0" 

CD 

CD 

03 

0 

0 

» 

5 

0 

0 

O 

0 

0 

qn 

I 

03 

03 

0 

0 

0 

CD 

E 

E 

CD 

0 

0 

GO 

D 

D 

CD  n  3  00 
if)  Zj  CO  CO 
LL  LL  CO  CO 
CO  CO  Q  Q 


"O  C 

C  o 

s  11 
Ills 

W  w  o 


C 
03 
CD 

^  £  £ 


0 

CD 

JO 

CD 

E? 

<0  3 

0) 

E  cr 

il 

O'  LL 


_ CD  CD 

Cg  0  O)  03 
(15  03  u  i- 

_Q  n  CD  CD 


0 

0 

03 

0 

0 

E 

be 

03 

03 

0 

_0 

0 

0 

0 

0 

cr 

0 

cr 

0 

1 

0 

1 

0 

03 

03 

1 

V 

0 

0 

CO 

£ 

E 

E 

CO 

0 

0 

0 

0 

LL 

LL 

LL 

O 

O 

LL 

0 

E 

0 

CO 

CD 

< 

z 

CO 

_J 

_i 

_ ! 

0 

LL 

LL 

CO 

CO 

0 

LL 

LL 

Q 

Q 

O' 

Output  File  -Sensitivity  Analysis 


o 

t  ■ 

o 

Cl 

T“ 

c 

C/) 

o 

o 

0 

o 

5= 

in 

c 

LU 

CO  q  3  CO 

to  J  w  to 
LL  LL  CO  CO 
CO  CO  D  Q 


•E  <D  W 
«  o)  a> 

£  te  fe 

O)  ©  © 

2'  O)  05 
CD  CD  CD 


0  ©  O 

i_  CD  CD 
V  CD  CD 

5  E  E 

O  (D  (C 
Cl  Q  Q 


.2  ^  CO  O  _j  < 
Q.  — I  CO  — I  -rj 
>  LL  LL  U_  CO  CO 

h  LL  LL  ii  Q  Q 


©C0Qr-T-0)^L0(DOC000 

COU)(DCDCOSOOOOO)C\J^ 

„(DO)CNjincoocoiocoNO)o 

Qcdcdssncocococococoo) 

<000000000000 

LUddddddddddod 


»2 


015  0.018  0.009  -0.014  0.004  -0.006  -0.002  0.003 

MinE  0.018  MinE50 -0.014  MinEIOQ  -O.OOfi  MinE250  0.003 
MaxEl0.033  I  MaxE50l0.005  I  MaxEIOQ  0.020  MaxE250  0.029 


Output  File  -  Sensitivity  Analysis 


Q) 

Q. 

O 

(rt 


c  =1 

O  S  O 
O  m 

T5  5= 

c  »-  c 

3  CD  - 

*  au 

Q-  C 
O  3  CM 

03  03 

®.EE 

<D  CO  CO 

!?  <D  CD 
CD  O)  03 
“  CD  CD 


$  $  CO  to 

50  i  B 

'v  'v  03  03 
CD  03  03  CD 

03  03  E  E 

S  £  CD  CD 

CO  CO  Q  Q 


CO 

CO 

LL. 

CO 


o  3  CO 
_J  CO  CO 
LL  CO  CO 
CO  Q  Q 


*0 
c 
3 

0  -o  ^ 

?  Ill 

to  *g  I  'n 

—  £  03 
CO  CO  03 


03 


CO 


'EE 

F  S? 
,  v>  £ 


03  03 
03  03 
CD  CD 


E  o- 

<D  0> 

2  & 

§  J 

a  Cl 

"S  11 


I  § 


03  03 
03  03 
CD  CD 

E  E 

CD  CD 

Q  O 


a  . 

>%  L 


CO  CD  H  < 

CO  _ I  — I  — ] 

u.  LL  CO  CO 
LL  LL  Q  Q 


E250 

03 

CM 

LO 

00 

c\i 

CM 

co 

CM 

0 

00 

c\i 

CO 

h- 

CM 

Is- 

CD 

CM 

CO 

cq 

CM 

CD 

CM 

LO 

CM 

LO 

CM 

CM 

CM 

r- 

00 

CM 

1 

03 

CM 

CM 

CD 

CM 

CM 

CM 

CM 

CO 

CM 

CM 

CO 

q 

CM 

CM 

O 

CM 

Is- 

03 

T~ 

CM 

q 

T— 

I 

CO 

00 

T— 

I 

co 

Is-, 

T— 

q 

CM 

0 

10 

CM 

CO 

00 

0 

0 

0 

0 

0 

d 

r- 

0 

0 

0 

0 

0 

d 

r- 

0 

0 

0 

0 

0 

d 

r- 

0 

0 

0 

0 

0 

d 

Is- 

O 

O 

a 

Is- 

O 

O 

O 

O 

O 

d 

h- 

0 

0 

0 

0 

0 

d 

Is- 

O 

O 

O 

O 

O 

d 

r- 

0 

0 

0 

0 

0 

d 

CD 

O 

O 

O 

O 

O 

d 

CD 

O 

O 

O 

O 

O 

d 

CO 

0 

0 

0 

0 

0 

d 

CD 

O 

O 

O 

O 

O 

d 

CD 

O 

O 

O 

O 

O 

d 

CD 

O 

O 

O 

O 

O 

d 

CD 

0 

0 

0 

0 

0 

d 

CD 

O 

O 

O 

O 

O 

d 

10 

0 

0 

0 

0 

0 

d 

LO 

O 

O 

O 

O 

O 

d 

in 

0 

0 

0 

0 

0 

d 

0.000005 

in 

0 

0 

0 

0 

0 

d 

in 

0 

0 

0 

0 

0 

d 

in 

0 

0 

0 

0 

0 

d 

LO 

O 

O 

O 

O 

O 

d 

O 

m 

CM 

LU 

c 

2 

0 

in 

CM 

LU 

X 

ro 

2 

0 

0 

LLJ 

| 

1 

I 

1 

1 

| 

1 

1 

1 

1 

1 

1 

1 

1 

| 

CM 

CM 

03 

0 

CM 

10 

0 

CM 

0 

CM 

h- 

q 

CO 

q 

I 

1 

0 

CO 

CO 

r- 

N- 

CM 

1 

CD 

O 

O 

O 

O 

O 

d 

CD 

O 

O 

O 

O 

O 

d 

CD 

O 

O 

O 

O 

O 

d 

CD 

O 

O 

O 

O 

O 

d 

CD 

0 

0 

0 

0 

0 

d 

CD 

O 

O 

O 

O 

O 

d 

CD 

O 

O 

O 

O 

O 

d 

CD 

O 

O 

O 

O 

O 

d 

CD 

O 

O 

O 

O 

O 

d 

CD 

O 

O 

O 

O 

O 

d 

CD 

O 

O 

O 

O 

O 

d 

CD 

O 

O 

O 

O 

O 

d 

10 

0 

0 

0 

0 

0 

d 

LO 

O 

O 

O 

O 

O 

d 

10 

0 

0 

0 

0 

0 

d 

LO 

O 

O 

O 

O 

O 

d 

LO 

O 

O 

O 

O 

O 

d 

LO 

O 

O 

O 

O 

O 

d 

LO 

O 

O 

O 

O 

O 

d 

LO 

0 

0 

0 

0 

0 

d 

LO 

0 

0 

0 

0 

0 

d 

3 

0 

0 

0 

0 

d 

3 

0 

0 

0 

0 

d 

3 

0 

0 

0 

0 

d 

3 

0 

0 

0 

0 

d 

0 

0 

LLJ 

C 

is 

O 

O 

LLJ 

X 

CO 

2 

O 

LO 

LLJ 

CO 

TT 

CM 

00 

CO 

CM 

LO 

CM 

03 

CM 

CM 

LO 

CO 

CM 

1 

CO 

CO 

c\i 

LO 

CO 

CM 

1 

1 

1 

CO 

CM 

| 

1 

03 

O 

CM 

m 

0 

CM 

0 

CM 

Is- 

q 

CO 

q 

10 

CO 

T— 

1 

03 

r- 

1 

1 

1 

I 

S50 

1 

LO 

O 

O 

O 

O 

O 

d 

to 

0 

0 

0 

0 

0 

d 

LO 

O 

O 

O 

O 

O 

d 

LO 

O 

O 

O 

O 

O 

d 

LO 

0 

0 

0 

0 

0 

d 

LO 

O 

O 

O 

O 

O 

d 

10 

0 

0 

0 

0 

0 

d 

10 

0 

0 

0 

0 

0 

d 

LO 

O 

O 

O 

O 

O 

d 

LO 

O 

O 

O 

O 

O 

d 

LO 

O 

O 

O 

O 

O 

d 

s 

0 

0 

0 

0 

d 

s 

0 

0 

0 

0 

d 

s 

0 

0 

0 

0 

d 

0 

0 

0 

0 

0 

d 

O 

O 

O 

O 

O 

d 

O 

O 

O 

O 

O 

d 

3 

O 

O 

O 

O 

d 

0 

0 

0 

0 

0 

d 

0 

0 

0 

0 

0 

d 

O 

O 

O 

O 

O 

d 

3 

0 

0 

0 

0 

d 

Tfr 

O 

O 

O 

O 

O 

d 

Tf 

O 

O 

O 

O 

O 

d 

CO 

0 

0 

0 

0 

0 

d 

0 

in 

LU 

c 

0 

in 

LU 

X 

(0 

LU 

1 

I 

1 

I 

1 

s 

CM 

r- 

Is- 

cm 

CM 

Is- 

CM 

03 

CO 

CM 

CD 

CD 

CM 

0 

CD 

CM 

S 

CM 

N- 

'=t 

CM 

CM 

CM 

h- 

00 

CM 

1 

1 

T— 

CM 

CM 

10 

CM 

0 

CM 

3 

CM 

0 

0 

CM 

a 

03 

00 

CO 

0 

co 

0 

CO 

CO 

0 

00 

1 

00 

0 

0 

0 

0 

0 

d 

00 

0 

0 

0 

0 

0 

d 

00 

0 

0 

0 

0 

0 

d 

00 

0 

0 

0 

0 

0 

d 

00 

0 

0 

0 

0 

0 

d 

00 

0 

0 

0 

0 

0 

d 

00 

O 

O 

O 

O 

O 

d 

r- 

0 

0 

0 

0 

0 

d 

Is- 

O 

O 

O 

O 

O 

d 

Is- 

O 

O 

O 

O 

O 

d 

Is- 

O 

O 

O 

O 

O 

d 

Is- 

O 

O 

O 

O 

O 

d 

Is- 

O 

O 

O 

O 

O 

d 

Is- 

O 

O 

O 

O 

O 

d 

CD 

O 

O 

O 

O 

O 

d 

CD 

O 

O 

O 

O 

O 

d 

CD 

O 

O 

O 

O 

O 

d 

CD 

O 

O 

O 

O 

O 

d 

CD 

0 

0 

0 

0 

0 

d 

CD 

O 

O 

O 

O 

O 

d 

LO 

0 

0 

0 

0 

0 

d 

in 

0 

0 

0 

0 

0 

d 

LO 

O 

O 

O 

O 

O 

d 

in 

0 

0 

0 

0 

0 

d 

in 

0 

0 

0 

0 

0 

d 

LU 

c 

2 

LL 

X 

TO 

^Parameter 

Value 

0 

0 

0 

LO 

0 

0 

0 

0 

0 

LO 

0 

0 

0 

CM 

0 

0 

LO 

CM 

0 

0 

0 

CO 

0 

0 

LO 

00 

0 

0 

q 

0 

0 

LO 

0 

0 

0 

LO 

0 

0 

LO 

LO 

0 

0 

0 

CD 

0 

0 

10 

CD 

0 

0 

0 

Is-' 

0 

0 

10 

Is-’ 

0 

0 

0 

00 

0 

0 

!  LO 

1  CO 

0 

0 

0 

03 

0 

0 

LO 

03 

0 

0 

0 

d 

0 

0 

10 

d 

0 

0 

q 

0 

0 

in 

r~ 

0 

0 

0 

CM 

0 

0 

in 

CM 

Parameter 

Value 

0 

0 

0 

LO 

CM 

O 

O 

LO 

CM 

0 

0 

q 

CM 

0 

0 

10 

CO 

CM 

O 

O 

O 

CO 

CM 

0 

0 

LO 

CM 

CM 

0 

0 

0 

cm' 

CM 

O 

O 

LO 

CM 

O 

O 

LO 

d 

CM 

O 

O 

O 

d 

CM 

O 

O 

LO 

03 

O 

O 

O 

03 

O 

O 

LO 

CO 

0 

0 

0 

06 

0 

0 

10 

h-' 

0 

0 

0 

Is-’ 

;  0 
0 

i  LO 
CD 

O 

O 

O 

CD 

O 

O 

LO 

in 

0 

0 

0 

LO 

0 

0 

in 

’'t 

0 

0 

q 

O 

O 

m 

CO 

T“ 

0 

0 

0 

CO 

0 

0 

in 

CM 

EAD250 

h- 

03 

CO 

O 

O 

d 

Is- 

LO 

CO 

O 

O 

O 

Is- 

CO 

O 

O 

d 

00 

r-. 

CM 

O 

O 

d 

03 

CO 

CM 

O 

O 

d 

CM 

CM 

0 

0 

d 

CO 

CO 

0 

0 

d 

CO 

CD 

8 

d 

03 

CO 

0 

0 

d 

CO 

r- 

0 

0 

0 

d 

0 

0 

0 

0 

0 

d 

O 

O 

O 

O 

O 

d 

O 

O 

O 

O 

O 

d 

0 

0 

0 

0 

0 

d 

0 

0 

0 

0 

0 

d 

O 

O 

O 

O 

O 

d 

O 

O 

O 

O 

O 

d 

O 

O 

1  O 
O 
!  O 

!  d 

O 

O 

O 

O 

O 

d 

0 

0 

0 

0 

0 

d 

O 

O 

O 

O 

O 

d 

0 

0 

0 

0 

0 

d 

0 

0 

0 

0 

0 

d 

O 

O 

O 

O 

O 

d 

0 

0 

0 

0 

0 

d 

0 

0 

0 

0 

0 

d 

EAD100 

<0 

CD 

03 

O 

O 

d 

00 

CD 

00 

O 

O 

d 

CD 

Is- 

r- 

0 

0 

d 

03 

0 

Is- 

O 

O 

d 

00 

CO 

CD 

0 

0 

d 

00 

CD 

LO 

O 

O 

d 

CD 

03 

s 

0 

d 

"3- 

CM 

0 

0 

d 

CO 

CO 

0 

0 

d 

CD 

CO 

CO 

0 

0 

d 

0 

03 

CM 

O 

O 

d 

LO 

rf 

CM 

O 

O 

d 

Is- 

0 

CM 

O 

O 

d 

03 

Is- 

T— 

O 

O 

d 

CM 

LO 

0 

0 

d 

0.00124 

0 

0 

0 

0 

0 

d 

0 

0 

0 

0 

0 

d 

0 

0 

0 

0 

0 

d 

0 

0 

0 

0 

0 

d 

0 

0 

0 

0 

0 

d 

0 

0 

0 

0 

0 

d 

0 

0 

0 

0 

0 

d 

0 

0 

0 

0 

0 

d 

0 

0 

0 

0 

0 

d 

0 

0 

0 

0 

0 

d 

EAD50 

h- 

CD 

CD 

O 

d 

00 

03 

CD 

O 

d 

Is- 

CM 

LO 

0 

d 

00 

00 

CO 

0 

d 

CD 

10 

CM 

0 

d 

CO 

CM 

O 

d 

0 

03 

03 

O 

O 

d 

CD 

CO 

00 

0 

0 

d 

is: 

r- 

0 

0 

d 

CM 

03 

CD 

O 

O 

d 

10 

5 

0 

0 

d 

00 

CO 

LO 

O 

O 

d 

00 

LO 

’<3- 

0 

0 

d 

CO 

00 

CO 

0 

0 

d 

0 

■t 

00 

0 

0 

d 

03 

CM 

O 

O 

d 

T— 

CM 

O 

O 

d 

r- 

03 

0 

0 

d 

r- 

CD 

O 

O 

d 

CO 

CO 

0 

0 

d 

CM 

0 

0 

0 

d 

0 

0 

0 

0 

0 

d 

0 

0 

0 

0 

0 

d 

0 

0 

0 

0 

0 

d 

0 

0 

0 

0 

0 

d 

0 

0 

0 

0 

0 

d 

EAD 

0 

03 

CD 

CD 

O 

d 

CM 

Is- 

CD 

O 

d 

T— 

Is- 

O 

CD 

O 

d 

CO 

r- 

co 

10 

0 

d 

Is- 

Is- 

CM 

LO 

0 

d 

00 

CO 

03 

"3- 

0 

d 

Is- 

03 

LO 

0 

d 

s 

CM 

tT 

O 

d 

CM 

03 

CO 

0 

d 

3 

LO 

CO 

0 

d 

CO 

T“ 

CO 

CO 

0 

d 

03 

CO 

0 

CO 

0 

d 

| 

5 

LO 

CM 

0 

d 

LO 

CM 

CO 

CM 

0 

d 

CO 

03 

O 

CM 

q 

d 

3 

03 

O 

d 

CM 

CM 

Is- 

O 

d 

CO 

LO 

LO 

O 

d 

10.01401 

LO 

10 

CM 

0 

d 

CO 

0 

d 

■<* 

03 

03 

O 

O 

d 

0 

03 

00 

0 

0 

d 

3 

r- 

0 

0 

d 

m 

0 

r- 

0 

0 

d 

103 


Output  File  -  Sensitivity  Analysis 


c 

o 

o 

c 

3 

4— 

(1)  O 

a.  > 

_0  _Q) 

c/>  <d 

1  ® 


03 
W  « 

5  g 

o  o 

V  V 
0  0 
O)  CD 

3  jS 
tt>  co 


?i 

II 

C 

0  — 
CL  "D 
Cl  C 
3  CM 


0  03 

E  E 

0  0 

0~  o' 

CD  03 
0  0 
-4— >  -4-4 

0  0 
I  I 
0  0 
CD  03 
0  0 

E  E 
0  0 
Q  Q 


CD  Q  3  CD 
co  Ij  co  CO 
LL  U_  CO  CO 
CO  CO  Q  Q 


> 

0 

"O 


■O  C 

c  o 
o  o 

E  ® 

1 1 1 

0  o  ^ 

.E  o'  a> 

0  CD  03 

E  S  « 

03  ®  ® 

b:  03  03 
0  0  0 

0  0 


CT 

O 


0  0 
03  03 
0  0 

E  E 
0  0 
E  Q  Q 


5 

o 


0  "  . 

.a  s  co  o  — j  < 

Q.  _J  CO  -J  -t 

>  LL  LL  LL  CO  CO 

1-llilu.QQ 


o 

LO 

CM 

LU 

l 

CO 

CM 

00 

CM 

h- 

03 

00 

00 

00 

CO 

■^r 

co 

CM 

O 

in 

CM 

CO 

1 

03 

CO 

o 

CO 

o 

00 

in 

o 

in 

in 

d 

o 

in 

CM 

LU 

c 

is 

O 

in 

CM 

LU 

X 

0 

S 

o 

o 

LU 

i 

CO 

in 

CM 

00 

00 

CM 

CO 

CM 

in 

o 

cm 

in 

o 

CM 

CO 

in 

CM 

O 

O 

CO 

1 

03 

CO 

o 

CO 

o 

00 

in 

o 

in 

in 

o 

o 

o 

LU 

c 

o 

o 

LU 

X 

0 

S 

§ 

LU 

h- 

03 

CM 

00 

Is- 

CM 

co 

in 

CM 

c\i 

T— 

c\i 

Is- 

03 

CM 

O 

in 

CO 

00 

CO 

d 

CO 

d 

00 

in 

o 

in 

in 

d 

o 

in 

LU 

c 

is 

O 

in 

LU 

X 

0 

LU 

CM 

CM 

in 

o 

CM 

co 

00 

Is- 

N; 

T“ 

Is- 

N- 

CM 

CM 

CO 

03 

CO 

d 

^r 

co 

d 

00 

in 

d 

in 

in 

d 

LU 

C 

is 

LU 

X 

0 

S 

3 

o 

E 

2  <d 

0  =l 

Q_  0 

<3  > 

o 

o 

o 

o 

o 

o 

CM 

d 

o 

CO 

o 

in 

co 

d 

L. 

0 

o 

1  ® 

2  3 

CD  CD 
CL  > 

CM 

CM 

d 

CM 

CO 

d 

CM 

d 

CM 

in 

o 

Is- 

m 

d 

o 

in 

CM 

Q 

< 

LU 

N- 

03 

CO 

o 

o 

o 

o 

o 

o 

o 

d 

00 

03 

CO 

o 

o 

d 

03 

03 

co 

o 

o 

d 

o 

o 

o 

o 

O 

O 

5 

< 

LU 

CO 

CD 

03 

O 

O 

o 

CO 

o 

o 

o 

d 

00 

03 

03 

O 

o 

d 

Is- 

03 

03 

O 

O 

d 

co 

03 

03 

O 

O 

d 

O 

in 

Q 

< 

LU 

h- 

co 

00 

o 

d 

CM 

o 

o 

CM 

o 

o 

^r 

o 

CM 

o 

d 

03 

03 

O 

o 

00 

00 

03 

o 

d 

Q 

< 

LU 

o 

03 

00 

CO 

o 

d 

00 

o 

00 

CO 

d 

Is- 

Is- 

03 

d 

Is- 

m 
co 
!  ^ 

!  ™ 
o 

in 

Is- 

CM 

CO 

CM 

o 

104 


Output  File  -  Sensitivity  Analysis 


0  O 

s-  > 
_0  0 

to  0 


i-i 

£  | 
u-  C 
0  ■“ 
CL  "D 
CL  C 
=5  CM 


|  E> 

OT  13 

§  § 

"f  V 

<D  0 
CD  CD 

0  3 

CO  CO 


0  CD 

E  E 

to  to 

CD  CD*' 
CD  CD 
CD  CD 
CO  to 

0  Cl) 
CD  CD 
CD  CD 

E  E 

CD  CD 

Q  O 


3  CO 
CO  CO 
CO  CO 
CO  CO  Q  Q 


CD  O 
CO 


> 

CD 

“O 

c  2 
0  10 

^  .c 

-  CO 
C  0 
to  -Q 
0  = 
-Q  0 

0  E 

CD  to 


“O  c 

c  o 
o  O 


1 1 
to  O 


to 


.E  0-  0 

to  CD  CD 

0i—i— 
_Q  0  0 


O)  ©  © 

E?  CD  CD 
_0  0  0 
to  to 

I  I 
0  0 
CD  CD 
0  0 

E  E 
0  0 
Q  O 


cr 

0 

o 


.52  S  w  CD  —I  < 

Q,  I  CO  — I  — J  rd 

>%  LL  LL  LL.  CO  CO 

h  LL  IL  LL  O  O 


E250 

0.0311 

o 

o 

CO 

o 

o 

0.0294 

0.0292 

0.0283  1 

0.0275  1 

0.0247 

0.0243  1 

0.0235  1 

0.0228 

CO 

CM 

CM 

O 

d 

0.0223  1 

0.0311 

S250 

o 

N- 

O 

O 

CO 

CD 

o 

d 

r- 

CD 

O 

d 

CD 

CD 

O 

d 

CD 

O 

d 

CM 

CD 

O 

o 

1 

CD 

LO 

O 

d 

LO 

LO 

O 

d 

CO 

LO 

o 

d 

CM 

LO 

o 

o 

LO 

p 

d 

MinE250 

MaxE250 

E100 

0.0297 

0.0289 

00 

CM 

o 

o 

0.0278 

0.0264 

0.0264 

1 

CD 

CO 

CM 

O 

o 

CO 

CM 

o 

d 

00 

CM 

CM 

o 

d 

CM 

CM 

O 

o 

0.0217 

0.0217 

h- 

CD 

(M 

o 

o 

O 

o 

CO 

CO 

LO 

o 

d 

CM 

LO 

O 

d 

LO 

o 

o 

o 

LO 

O 

d 

p 

d 

o 

o 

1 

CM 

o 

o 

CM 

o 

d 

T— 

o 

d 

o 

'M- 

o 

o 

CD 

CO 

o 

d 

o 

o 

m 

c 

o 

o 

LU 

X 

0 

:> 

0S3 

0.0236 

92Z00 

0.0224 

0.0223 

T— 

CM 

O 

o 

0.0216 

0.0190 

0.0193 

00 

oo 

o 

d 

B 

o 

CO 

o 

d 

■5 

CD 

CO 

CM 

O 

o 

■ 

o 

CO 

o 

d 

CD 

CM 

O 

d 

00 

CM 

O 

o 

CO 

CM 

O 

d 

CM 

o 

o 

N- 

CM 

O 

d 

CM 

O 

d 

CM 

O 

o 

CM 

o 

d 

CO 

CM 

O 

o 

CO 

CM 

o 

d 

MinE50 

o 

LO 

LU 

X 

0 

LU 

S6000 

0.0095 

0.0099 

CM 

o 

o 

d 

66000 

0.0105 

1 

Z6000 

o 

p 

d 

o 

o 

o 

d 

0.0100 

o 

p 

d 

LO 

CD 

O 

O 

o 

LO 

o 

O 

o 

CO 

O 

CO 

o 

d 

o 

CO 

o 

d 

CO 

o 

o 

CM 

CO 

o 

d 

CO 

o 

d 

CO 

CO 

o 

o 

1 

O 

CO 

o 

o 

CO 

o 

o 

CO 

o 

d 

T 

CO 

o 

o 

CM 

CO 

o 

d 

LU 

c 

LU 

X 

0 

AParameter 

Value 

o 

CD 

d 

o 

LO 

d 

o 

d 

o 

CO 

d 

o 

CM 

d 

o 

d 

o 

o 

d 

o 

o 

i 

o 

CM 

d 

1 

o 

o 

i 

o 

LO 

d 

i 

Parameter 

Value 

p| 

I 

o 

o 

d 

o 

d 

o 

CM 

d 

o 

CO 

o 

o 

d 

EAD250 

0.06103 

CD 

CO 

CD 

CD 

O 

d 

LO 

O 

O 

o 

0.07506 

0.07876 

0.08241 

0.08539 

0.08795 

0.09014 

0.09200 

0.09358 

0.09487 

o 

o 

Q 

< 

LU 

1 

S022L0 

0.12501 

0.12759 

0.12974 

0.13153 

0.13299 

0.13418 

T~ 

LO 

CO 

d 

0.13582 

CM 

CO 

CD 

CO 

d 

0.13659 

EAD50 

0.18510 

0.18575 

0.18616 

CO 

CM 

CD 

00 

d 

0.18605 

0.18572 

0.18520 

N- 

LO 

00 

o 

oo 

CO 

00 

o 

0.18294 

1 

0.18078 

EAD 

0.32974 

10.32679 

0.32423 

0.32142 

o 

00 

00 

d 

0.31516 

O 

CD 

T— 

CO 

d 

CD 

00 

00 

o 

co 

o 

CO 

CD 

LO 

O 

CO 

d 

0.30251 

1 0.29940 

0.29608 

105 


Output  File  -  Sensitivity  Analysis 


CD  Q  Z>  CD 
05  Zj  CO  CO 
Li_  LL  CO  CO 
CO  CO  Q  Q 


"O  g 
C  o 

s  11 

?  1 1£ 


CO  CD  CD 
CD  03  *-  *- 

-Q  _q  CO  CO 


(P  CD  (C 

-  W  (S 


rr  rT  1  * 

0)  0  ©  © 

C.  C.  CD  CD 

O 

V  V  CD  CD  CQ 

to 

£  £  E  E  co 

o 

O  O  CO  03  LL 

<c 

LL  LL  Q  Q  CO 

UJ 

CD 

E 

CO 

CO  e)  — 1  <  z 

co  _j  c 

Q 

LL  LL  CO  CO  3 

< 

LL  LL  Q  Q  C£ 

LLI 

16 


Output  File  -  Sensitivity  Analysis 


m  O  3  CD 
CO  _J  CO  CO 
UL  LL  CO  CO 
CO  CO  O  O 


5  -S I 

■o  5  fc  C 

2  8  I  - 

<0  «  o  “ 

.?=  ci 


CO  o  -I  < 
CO  Zj  — 1  — 1 
LL  il  CO  CO 
LL  LL  Q  Q 


17 


Output  File  -  Sensitivity  Analysis 


CO  o  ID  CD 
C/3  _j  C/D  C/D 
LL  LL  C/D  C/D 
C/D  C/D  D  Q 


"D  £ 
C  o 

o  "S 

>  u. 

5  8™ 

(0  O  ^ 

.E  0"  CD 
(/)  CD  O) 

5  0  0 
§>  ®  o* 

uD)  D) 
0  0  0 
-  00  « 

®  i  i 

1-  03  O) 

V  0  0 
|  E  E 
o  0  0 
LL  O  O 


o  ^  C/D 
a  j  c/D 

>1  LL  U- 


0.  C/D  C/D 

ll  a  a 


Output  File  -  Sensitivity  Analysis 


L_ 

Q) 

Q 

E 

ffi  o  2  “ 

TO  a) 

CO  _ I  C/5  CO 

CO  3 

u_  ul  co  co 

Q_  CO 

CO  CO  Q  Q 

<  > 

■o  5  2  “ 
to  m  |  » 

■—  .E  CD  < D 

<2  co  o)  05 

(0  OJ  1-  L. 

-O  .q  to  CO 


CO 

— 

CD 

n 

CO 

O) 

Q) 

CD 

CD 

CD 

E 

CO 

u_ 

JO 

CD 

ro 

CD 

_C0 

v_ 

00 

« 

CO 

cr 

cr 

a> 

CD 

i 

CD 

CD 

_J 

CD 

L. 

(D 

L_ 

tx 

M— 

1 

1 

CO 

CO 

CD 

£ 

E 

E 

4= 

o 

o 

CO 

CO 

1 

1 

LL 

LL 

O 

Q 

C /)(!)  J< 
CO  — I  — I  rd 
1 1  1 1  CO  CO 

LL  LL  O  Q 


>9 


Output  File  -  Sensitivity  Analysis 


DQ  Q  3  CD 
(f)  y  CO  CO 
LL  U.  CO  CO 
W  CO  Q  Q 


S  -i  J 

E  cr  V 
ffi  ®  5 
•y  *+=  o 


w  _o  ™ 

.E  <D  CD 
w  O)  D) 

5  te  fe 

a  <d  © 

CJ)  CD 
_CD  CD  CD 

-  00  00 
rj  •  1 

o  © 
£  E  E 

O  05  CD 
LL  Q  D 


«  2  co  0  —I  < 

Q. _ I  CO  — I  rd 

>  LL  LL  LL  JO  CO 
hu.iLU.QQ 


Output  File  -  Sensitivity  Analysis 


c 

o 

c 

o 

u 

o 

0 

o 

-Q 

c 

c 

13 

0 

■  — 

M— 

CL 

T3 

bf— 

CL 

C 

0 

o 

3 

CM 

CL 

> 

„ 

O 

CD 

0 

0 

CO 

0 

E 

E 

— 

<D 

CO 

co 

03 

CD 

p. 

E 

i_ 

0 

0 

£/) 

03 

CD 

CD 

_ 

0 

0 

£ 

5 

*0 

0 

o 

o 

0 

0 

| 

i 

CD 

CD 

0 

0 

0 

0 

CD 

CD 

E 

E 

0 

0 
-♦— > 

0 

0 

CO 

CO 

Q 

Q 

CQ  Q  Z>  CO 
(O  _J  CO  CO 
LL  U_  CO  CO 
CO  CO  Q  Q 


> 

0 

■O 

5 

“O 

0 

j*: 

CO 

0 

c 

C 

0 

0 

0 

0 

-Q 

-Q 

— 

0 

CD 

E 

0 

i_ 

0 

CT 

CT 

0 

E 

B: 

o 

o 

LL. 

LL 

E 

E 

CO 

O 

0 

0 

CO 

< 

Q 

Q 

Q 

LU 

0 

E 

0 

_j 

< 

2 

_i 

c 

Q 

CO 

CO 

3 

< 

Q 

Q 

DC 

LU 

Output  File  -  Sensitivity  Analysis 


C 

o 

o 

c 

3 

H— 

(D  o 

9-  > 

_o  _0) 

CO  <D 
==  CD 
| 

0  « 

3:  3 
o  o 
5=  5= 
CD  0 
03  CD 
03  CO 

CO  CO 


Y  | 

1  | 
s-  C 
(D  — 
Q.  "O 
Q.  c 
3  CNJ 


03  03 

E  E 

CO  CO 

CD  CD~ 
CD  CD 
03  03 

t 0  to 

CD  CD 
CD  O) 
03  CO 

E  E 

03  CO 

Q  0 


CO  Q  Z)  CO 
CO  -J  CO  CO 
u_  LL  CO  CO 
CO  CO  Q  Q 


> 

0 

■O 

"O 

0 


c 

*0 

0 

_D 

0 

03 

l_ 

//I  ro 

0  _j 

g  - 

E  o" 

z  I 

il 

O'  LL 


■O  c 

c  .2 

11 

ill 

0  O  ic 


■E  .2  <D  0~ 

52  0  03  03 

0  CO  21  i- 

n  _q  0  0 


=  0 


0  ^  0  CD 

C  E?  03  03 
_0  0  0 
-  0  0 

S"  S"  ^  ^ 

V  -*r  TO  « 

5  5  E  E 

O  O  0  0 
LL  LL  Q  Q 


a  -i 

>  u_ 


co  O  -!  < 

CO  _l  — 1  — 1 

ul  LL  CO  CO 
LL  U_  O  Q 


2 


Output  File  -  Sensitivity  Analysis 


c 

.2 

C 

3 

o 

0 

0 

m— 

U— * 

o 

-Q 

c 

c 

3 

0 

*“ 

M— 

CL  "O 

M — 

CL 

c 

0 

O 

=3 

CM 

CL 

> 

O 

0 

0 

0 

0 

E 

0 

E 

*0 

E 

<D 

0 

0 

CO 

L_ 

0“ 

0 

0 

0 

CO 

CO 

__ 

0 

TO 

3? 

0 

0 

o 

o 

0 

0 

V 

*r 

CO 

CO 

0 

0 

0 

0 

O) 

0 

CO 

0 

E 

0 

E 

0 

00 

CO 

O 

Q 

CQ  O  3  CD 

(/)  J  W  (/) 

LL  LL  CO  CO 
CO  CO  o  o 


s  § 
s  11 

^  5  fe  C 
OO  M  I  “ 


•,t:  •*=  ® 


b  D  i  ro  CD 


c 

‘co 

CD 

-Q 


a  '8  *» 


0 

CO 


CO  S)  0  0 
F  ?  CD  CD 

w  J5  2  2 

0  0 


0 


©  0  <i)  o 


TO  g 

Z  V 


3 


11 

LL  LL  O  O 


O)  CO 
CO  CO 

E  E 

CO  0 


m  II 


Q.  _J 
>  U- 


co  o  j  < 

CO  _J  — l  — J 
LL  LL  CO  CO 
LL  LL  Q  Q 


3 


Output  File  -  Sensitivity  Analysis 


CQ  n  3  CQ 
(f)  Zj  CD  CO 
U_  LL  0)  CO 

0)0)0  o 


c 

2 

0 

0 

5 

(/) 

03 

0) 

« 

0 

o 

c 

c 

0" 

0~ 

CO 

'c/3 

03 

03 

c 

03 

i— 

CO 

_Q 

_Q 

0 

_0 

03 

— 

0 

- 

-Q 

0 

O) 

03 

E 

to 

E5 

JO 

0 

O) 

0 

0 

03 

0 

i_ 

„ 

0 

ts 

03 

D" 

O’ 

i 

0 

03 

1 

0 

03 

— 1 

0 

i— 

£ 

B  2  5 

7  ^  o 


o  ^  co 

Q.  _l  O) 
>s  U-  LL 
I—  LL  LL 


I  0  CO 

5  E  E 

O  CO  CO 

EOQ 


LL  CO  O) 

LL  O  O 


Output  File  -  Sensitivity  Analysis 


c 

o 

o 

c 

13 

**- 

0  O 

O  _Q> 

0  03 

^  © 


03 


?  -I 

II 

c 

0)  -- 
CL  -O 
Q-  C 
D  CM 


03  03 

E  E 

0  CO 


E 

Vr  rrt 


§  I 

"F  **? 

CD  CD 
CD  CD 
03  03 

CO  CO 


P>  m 


0  0 
CD  CD 

0  0 
-*— ■  -4— « 

0  0 
i  i 
0  0 
CD  CD 
0  0 
E  E 
0  0 
Q  O 


3  CO 
CO  CO 
CO  CO 
CO  CO  Q  Q 


CO  O 
CO 


T3  £ 

c  .2 
O  O 
2  © 


^  5  fc  C 
u  0  > 

»  0  o 


0 


0 

0 

_Q 

0 

CD 

l_ 

M  03 

CO  _ I 

<D  -1 

E  cr 

0  © 

Z 


•§  .E  <d  0“ 

£  0  CD  CD 

0  (0  L  t- 

JQ  _Q  jo  _0 

i  ra  (D  fl) 

C  E?  CD  03 

0  _E  j2  iS 
-00 

S"  S"  ^  ^ 

V  V  0  0 
|  |  E  E 

jo  _o  0  0 

LL  LL  Q  D 


c  & 

3  O 
O'  LL 

"0  11 

O  ^  (O  0  J  < 

5.  =i  co  Ij  -j  -J 

>»  LL  LL  LL  00  CO 
h  L  LL  LL  O  O 


E250 

ozs 

5.30 

N- 

00 

4.36  I 

h*- 

00 

co 

2.89 

2.55 

2.25 

ID 

CM 

CM 

0 

h-*- 

LO 

S250 

0.029 

0.027 

0.025 

2200 

1 

E 

SLO'O 

0.013 

O 

d 

O 

LO 

CM 

LU 

c 

0 

LO 

CM 

LU 

X 

0 

E100 

LO 

CM 

CD 

co 

LO 

CO 

LO 

00 

r^- 

1 

03 

CD 

CO 

CD 

CO 

03 

h- 

CM 

LO 

CM 

LO 

CM 

LO 

CM 

CD 

S100 

6200 

0.027 

0.025 

0.022 

0.017 

0.015 

CO 

0 

d 

O 

d 

O 

O 

E 

c 

O 

O 

LU 

X 

0 

E50 

CD 

CO 

CO 

CO 

CD 

h- 

CNJ 

CD 

0 

CD 
i 0 

1 

LO 

CM 

03 

LO 

CO 

CM 

CO 

h- 

CD 

CM 

h- 

CD 

CM 

CD 

CO 

S50 

0.029 

0.027 

0.025 

0.022 

1 

0.017 

O 

d 

CM 

O 

d 

T — 

O 

d 

O 

LO 

LU 

c 

0 

LO 

LU 

X 

0 

LL! 

CO 

LO 

03 

03 

00 

LO 

B 

1 

00 

T— 

CO 

CM 

h- 

CM 

0 

CM 

CM 

CM 

CM 

T — 

c\i 

h- 

CO 

LO 

CO 

0.029 

0.027 

0.025 

0.022 

1 

0.017 

0.015 

0.013 

O 

d 

LU 

c 

LU 

X 

0 

^Parameter 

Value 

O 

O 

O 

O 

O 

O 

CO 

2.000 

O 

O 

O 

O 

O 

O 

d 

O 

O 

q 

V— 

-2.000 

O 

O 

O 

CO 

t 

0 

0 

0 

Parameter 

Value 

O 

O 

CD 

CO 

0096 

O 

O 

CD 

d 

11.600 

12.600 

13.600 

O 

O 

CD 

0 

0 

CD 

LO 

0 

0 

CD 

CD 

EAD250 

0.00397 

0.00397 

0.00397 

0.00397 

0.00397 

CD 

03 

CO 

O 

O 

d 

CD 

O) 

CO 

0 

0 

d 

0.00396 

0.00395 

EAD100 

0.00977 

0.00975 

0.00973 

0.00971 

0.00968 

0.00964 

0.00960 

0.00955 

0.00940 

EAD50 

T- 

03 

O 

d 

S06100 

0.01897 

00 

CO 

00 

0 

d 

0.01866 

0.01819 

0.01746 

0.01645 

0.01487 

O 

< 

LU 

1^ 

00 

d 

0.14952 

0.09060 

CM 

CO 

00 

CD 

q 

d 

0.05109 

O 

CO 

00 

CO 

0 

d 

0.02927 

10.02240 

115 


Output  File  -  Sensitivity  Analysis 


S  J 


(D  O 

g-  > 

-9  _Q) 

0  0 


1  ® 
1-  c 
0  — 
CL  -O 
CL  C 
=5  C\1 


—  0 
| 


CO  CO 

E  E 

CO  CO 


CO 


0  0 
CO  CO 
0  0 


I  I 


CO  0 


0  0 
CO  CO 
0  0 


CO  CO 


0  0 
CO  CO 
0  0 
E  E 
0  0 
Q  O 


OQ  Q  3  CD 
CO  _J  CO  CO 
LL  LL  CO  CO 
CO  CO  O  O 


-O  c 

C  O 

o  o 

_Q  0 


c 

0 

0 


55  ■£ 
§2 


c  c 


_ 0  0 

£  CO  CO  CO 


0 

0 

-Q 


0  0  _  _ 
_Q  _o  _co  _ro 


0  o)  ©  © 

c  co  co 

^000 
0  CO 


0 

0 

E  cr 

ro  © 


0  0  C)  JP 


CO  CO 
0  0 

5  £  E  E 


§  1 

O'  U- 


O  O  0  0 
E  Q  Q 


o 


a  -i 

>  u- 

f—  u_ 


co  CD  H  < 
co  _J  — I  — I 

ul  lZ  CO  CO 

LL  LL  O  O 


6 


Output  File  -  Sensitivity  Analysis 


c 

o 

c 

3 

O 

o 

£5 

c 

u_ 

3 

0 

M— 

CL 

M— 

CL 

0 

O 

3 

CL 

> 

,2 

0 

0 

0 

0 

E 

:= 

0 

0 

0 

E 

05 

i_ 

0^ 

0 

0 

05 

““ 

0 

5 

0 

o 

o 

0 

V 

K 

05 

0 

0 

0 

05 

05 

E 

0 

0 

0 

CO 

35 

Q 

§  § 


O  o 
0  0 


c  C 


Y  1 

CM  CM 


I  E> 


0  0 
05  05 
CO  CO 


CO  CO 


0  0 
05  05 
CO  CO 

E  E 

CO  CO 

Q  Q 


CD  O  3  CO  CD 

CD  -J  CD  CD  _l 

LL  LL  CD  CD  CD 

CD  CD  O  Q  O 


> 

0 

~o 


■O 

c 

3 

o 


c 

o 


C 

0 

0 

0 

0 

C 

'0 

C 

*0 

C 

0 

X 

0 

X) 

*0 

0 

XI 

0 

E 

0 

05 

0 

0 

_0 

05 

0 

CJ 

O' 

0 

0 

0 

i_ 

0 

i— 

E 

0 

CT 

0 

i_ 

1 

5: 

0 

1 

O 

z 

V 

LL 

LL 

c 

3 

O 

CE 

LL 

0 

II 

0 

CD 

CD 

Q.  —I 

CD 

_l 

>  Li- 

LL 

LL 

LL 

LL 

LL 

o  s 


.Q  ® 


0 


.=  00 


0  0 
05  05 


05  05 
CO  0 


0  0 


llQQ 


7 


Output  File  -  Sensitivity  to  Uncertainty 


COO)^OCOCD(DLOLOCO^ 
oC0t~t—  t-OOOOOOO 
oOOOOOOOOOOO 

^qoqqqqqpqqo 

C/Dddddddddddd 


OLOCMCNt- 

lo  o  o  o  o 

LU  dodo 


h-  oo  in  co  O)  s  co 
^OOtDincOCMCOT-T-.CM 
T-OOOOOOOO^O 
oOOOOOOOOOOO 

Lopqqqqqqqqujq 

C/DOOOOOOOOOt-O 


(Noosinomco^co 

rOOOO^-OCMOOCOO 

OCDLO^-'sfCOCNJCMCNJ 

ooooooooo 


O)COCDN^-SIOCOt-O)C0 

OLOCOCvJCMt-x-t-t-OO 

T-OOOOOOOOOO 

oooooooqqqq 

coddddddddddd 


-=oooooooooo_ 

tDOJCONCDin^COCMr-OO 


—  O  O  O 

CC  CCoOOOOOOOOOLOO 
Q.>r-C\ICO^in0SOOO)T-T-C\! 


<0  O  CNJ  CN  Lfi  T-O)  LO^-t-O)  CD 
IjjCOOOCOlfi^^COCOfOfOCM  CM 
<1  CO  S  N  N  N  NN  MMS  s  N 


<Ot-OCOO)CDLOCOCNOOI^ 
LUCOOQOOSNSNNSCDCq 
<]  NNtD0CDCdcDcd(dcdc£)(£i 


<CNO)LOCNJOOOOOCDOUOLO^r 
^OOOOOOONNSNSSNN 
<  dduoLOLOinLOLodLOLdd 


QSCD^-fO^-^OinCMO^O 

<r^coooT-oa)qqqis-N:r--; 

lUoiodcdcdoossNSSNS 


;30  lMaxE5Q IQ.Q513  iMaxEIOQ  10.0994  lMaxE250 


Output  File  -  Sensitivity  to  Uncertainty 
Flow  Frequency  Function 


CD 

C 

0 


o 

o 

9? 

a> 

CD 


TO  — 


W  S' 

E  % 

TO  > 


> 

0 

■Q 

4 1 

if 
0  0 
E*  « 

TO  CD 
—  TO 

S 

O 
5= 


C 

3 

on 

TO 

.2  2 
a  -i 
>  ll 
H  ll 


c/3 
0 

.  CD 
0  TO 
CD  £ 

jS  TO 

co  a 


LL  CO 
CO  Q 


0 

E 

TO 


C 

3 

cr 


X“ 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

■ 

00 

O 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

LO 

CM 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

It 

d 

r— 

LU 

c\i 

03 

r- 

LO 

CO 

CO 

CM 

CM 

c\i 

T- 

■ 

E 

o 

CO 

CO 

CO 

CO 

CO 

CO 

o 

LO 

o 

LO 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

CM 

CM 

LU 

X 

O 

lO 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

UJ 

LU 

LU 

LU 

1 

LU 

CM 

CO 

03 

o 

CD 

CO 

X— 

CM 

LO 

o 

q 

TO 

CO 

LO 

CM 

CM 

X — 

03 

CD 

to 

h- 

03 

TT 

O 

00 

o 

00 

CM 

CM 

r- 

o 

T— 

CO 

'M- 

CD 

00 

CO 

03 

00 

CO 

o 

r- 

X— 

o 

X — 

CD 

00 

CM 

CM 

X— 

X— 

X— 

X — 

o 

o 

X— 

X“ 

T— 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

x— 

LU 

CO 

o 

o 

d 

o 

o 

o 

o 

o 

d 

o 

d 

d 

00 

CM 

o 

00 

CD 

LO 

LO 

CO 

CM 

1 

o 

o 

o 

o 

o 

CO 

x— 

x— 

X— 

o 

o 

o 

o 

o 

o 

o 

UJ 

LU 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

X 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

TO 

CO 

o 

d 

d 

o 

o 

o 

d 

d 

d 

d 

o 

CD 

LO 

00 

00 

CO 

CO 

03 

00 

03 

00 

o 

1 

o 

CO 

CD 

CO 

CM 

00 

CM 

03 

o 

00 

Tt 

o 

o 

CD 

o 

LO 

00 

CM 

X— 

x — 

o 

T— 

o 

o 

o 

■ 

o 

LO 

LO 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

■ 

o 

o 

LU 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

■ 

d 

d 

CM 

00 

CO 

CO 

03 

CO 

LO 

LO 

■ 

r — j 

o 

CO 

r- 

LO 

CO 

CM 

CM 

CM 

o 

o 

■ 

LO 

LU 

LO 

v— 

o 

o 

o 

o 

o 

o 

o 

1 

LU 

o 

UJ 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

X 

lO 

o 

o 

o 

o 

o 

o 

o 

o 

LU 

o 

TO 

CO 

d 

d 

d 

o 

d 

d 

d 

d 

o 

CO 

LO 

CO 

h- 

1^ 

o 

03 

00 

00 

CM 

CM 

CO 

00 

LO 

03 

CM 

CM 

Is- 

o 

r- 

00 

00 

CM 

LO 

o 

r- 

CD 

LO 

CO 

CO 

CM 

X — 

X — 

CM 

CO 

T— 

o 

o 

o 

o 

o 

o 

o 

o 

o 

CO 

LU 

o 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

CM 

LO 

CD 

X“ 

Is- 

LO 

CO 

o 

03 

CD 

o 

LO 

CO 

CM 

CM 

T“ 

o 

o 

UJ 

LU 

X— 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

X 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

c 

TO 

CO 

o 

d 

d 

d 

o 

d 

d 

o 

d 

d 

o 

i_ 

0 

0 

E 

<0 

k. 

0 

TO 

13 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

CL 

TO 

03 

00 

r- 

CD 

LO 

'M- 

CO 

CM 

X — 

o 

o 

< 

> 

V“ 

-H 

-H 

-H 

LO 

o 

k_ 

-2 

0 

E 

0 

5 

3 

o 

o 

o 

TO 

TO 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

LO 

o 

Q_ 

J> 

x — 

CM 

CO 

LO 

CD 

r- 

00 

03 

CM 

O 

LO 

CM 

a 

< 

LO 

r^ 

03 

o 

o 

r- 

M- 

CM 

o 

Is- 

LO 

LU 

B 

CD 

CO 

CO 

CM 

CM 

CM 

CM 

x— 

T“* 

<_ 

CD 

cb 

CD 

CD 

CD 

d 

CD 

CD 

CD 

CD 

CD 

O 

o 

Q 

< 

LO 

00 

Is- 

CM 

00 

LO 

00 

CM 

03 

r- 

CD 

LU 

00 

Is- 

CD 

CD 

CD 

CD 

E 

LO 

LO 

LO 

<_ 

E 

LO 

LO 

LO 

LO 

LO 

LO 

LO 

is 

LO 

LO 

LO 

O 

LO 

a 

< 

o 

03 

Tt 

CM 

o 

03 

00 

00 

CD 

CD 

LO 

LO 

UJ 

03 

Is- 

h- 

Is- 

r- 

CD 

CD 

CD 

CD 

CD 

1  CD 

CD 

<_ 

X5t 

xfr 

Tf 

a 

CM 

CO 

00 

o 

00 

00 

Tl- 

o 

00 

CM 

03 

< 

CD 

LO 

B 

o 

03 

00 

r- 

Is- 

r- 

CD 

CD 

LO 

LU 

CO 

1^ 

is 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

119 


Output  File  -  Sensitivity  to  Uncertainty 


Output  File  -  Sensitivity  to  Uncertainty 
Stage  Flow  Function 


CD 

C 

CD 


0 

"O 

>  2 

-g  ^ 


o 

o 

£ 

CD~ 

CD 


CT 

CD 


CD  E 

CD  W 

E3  CD 
CD  CD 
—  CD 

o 


ro  £ 

Z  o 

C  Ll_ 

3 


CD 

*v  2^ 

CD  CD 
CD  £ 
2  <6 
CO  Q 


<o 

2  z 

Q.  -J 

>  LL 

H  ll 


3  to 

LL  CO 
CO  Q 


CO 

CO 

LL 

CO 

CD 

E 

CD 

z 

c 

3 


Output  File  •  Sensitivity  to  Uncertainty 


c 

o 

43 

O 

c 

3 

Li- 

CU 

O) 

<0 

CO 

0) 

O) 

(0 

E 

to 

O 


05 

C 

05  > 

“  05 

2  -O 
O  -O 

CJ  1 

05  « 

i— 

-  05 
05  O) 

E?  0 

_0  “ 

» si 


0 

"O 

2 

GO 

0 

E 

CO 

CD 

05 

CO 


V  0 

o  jo 

U-  CO 


0) 

E 

to 

z 

c 

3 


to 

.2  z: 

Q.  — I  — 1 
>11-4- 
h-  U_  CO 


CO 


0 

05 

CO 

E 

0 

Q 


CO 

CO 

CO 

Q 


CO 

CO 

CO 

o 

0 

E 

CO 


c 

3 

DC 


Output  File  -  Sensitivity  to  Uncertainty 
Damage  Stage  Function 


CD 

C 

CD 


> 

CD 

*u 

>  X) 
%  * 


o  • 

o 

a) 

l_ 

CD 

O) 


2  2>< 

il 


c  “■ 

3 
O' 

I  z  CO  W 


Su. 


CO 
CO 
CO  o 


Output  File  -  Sensitivity  to  Uncertainty 


C 

o 


o 

c 

3 


5 

o 


0) 

o> 

<0 


CO 


D) 

c 

©  > 
“  CD 

V  -o 


> 

CD 

"D 

"D 

0 


■O 

3 


CO 


0 

>  03 

■g  O 

"O 


O  "O 

0  w 


0  0 
~  3 


0  P° 


-  0 
0  0) 
CD 

i_  0 
_0  ~ 

CT  % 
©  £ 


E  © 
w  9- 


_  o 
0  £ 


CO 

0 


w  £ 
0  .9 
O)  V 
0  0 


c 

3 

a: 


D) 

E 

CD 

0 

0 

0 

CO 

o 

CO 

0 

Cl 

O 

CO 

CO 

£ 

_l 

CO 

CO 

u_ 

CO 

LL 

CO 

Q 

co 

Output  File  -  Sensitivity  to  Uncertainty 
Stage  Flow  Function 


> 

0 

“O 


2 

to 


T3 

3 


CO 


<D 

(/) 

CD 

o 


>  § 
■S  £ 
>  2  & 

ip 

■“  2  c 


(/) 


0 

CD 

CT> 

i— 

t_ 

0 

JD 

cr 

O 

0 

L_ 

M— 

V 

0 

£ 

CD 

o 

0 

LL 

CO 

w  E  o 
0  «  E 

CD  CD 


D) 


3 


CD 

.y  ^ 

Q.  —I 

>  ll 
H  LL 


C 

o 

E 

CO  W  CD 
-J  W  X 
LL  CO  LL 
CO  Cl  CO 


c 

o 

E 


CD 

X 


CO 


0 

E 

CD 

2 


c 

3 

X 


25 


Output  File  -  Sensitivity  to  Uncertainty 


05 

c 

0 


O 

o 

0 


(/> 

0 

E 

0 

Z 

c 

3 

a: 

"ro 

o 


CD 

c 
0  > 
__  0 
■O  T3 

O  -Q 
8  « 
^  0" 

0  CD 

E>  fe 

0  — 

y| 

th  ^ 

V  0 

o  5 

LL  (/) 


> 

0 

~o 

T3 

«*-* 

0 


> 

X5 

3 

<75 

0 

0 

0 

a 


0 

0 

CL 

CL 


CO 

Q_  — I  — 1 

>  u~  4- 

I-  U.  CO 


0  O 

E  C 

0  O 

0~  ^ 
CD  0 

to  X 
0  - 

0  CD 
CD  £ 

0  V 
E  3S 
0  o 

Q  LL 

c 

o 

E 

CO  X 

co  £ 

O  LL 


Run  Name  FFHarmon 

O 

LO 

CM 

LU 

05 

o 

d 

LO 

Is- 

q 

d 

CM 

CD 

o 

d 

0.023 

LO 

q 

d 

00 

o 

d 

0.012 

o 

o 

d 

r— 

o 

o 

d 

h- 

o 

q 

d 

00 

o 

o 

d 

1 

O 

LO 

CM 

CO 

CM 

CM 

d 

CO 

iq 

CM 

LO 

d 

CM 

00 

Is- 

00 

q 

LO 

Is- 

CM 

x— 

q 

CD 

CM 

CO 

CD 

CM 

05 

00 

CM 

00 

CD 

c\i 

1 

o 

o 

LU 

Is- 

LO 

T— 

o 

CM 

LO 

O 

d 

00 

o 

d 

CM 

o 

d 

o 

CM 

o 

d 

CD 

X“ 

o 

d 

o 

d 

05 

o 

o 

d 

05 

o 

o 

o 

r- 

o 

o 

d 

00 

o 

o 

d 

1 

o 

o 

CO 

05 

d 

^r 

q 

CD 

00 

00 

d 

T“ 

LO 

CD 

CM 

CD 

CD 

o 

LO 

00 

CD 

05 

CM 

CD 

05 

CM 

CO 

CM 

LO 

CD 

CM 

1 

o 

LO 

LU 

00 

05 

q 

d 

o 

O) 

CO 

o 

o 

T“ 

CD 

CM 

o 

d 

CD 

05 

q 

d 

^r 

T— 

o 

o 

CM 

q 

d 

00 

05 

o 

q 

d 

CD 

CD 

o 

o 

d 

r- 

o 

o 

d 

CM 

oo 

o 

o 

o 

oo 

o 

o 

o 

d 

1 

O 

LO 

CO 

Is- 

00 

o 

CO 

CM 

CM 

Is- 

i^ 

X~ 

CD 

05 

00 

Is-’ 

o 

CM 

O) 

LO 

CD 

CD 

00 

oo 

N; 

00 

o 

CD 

05 

CM 

00 

o 

o 

CM 

o 

CM 

Is- 

LO 

05 

d 

CM 

o 

d 

1 

LU 

LO 

CM 

d 

CD 

LO 

o 

d 

05 

CO 

o 

d 

o 

CO 

o 

d 

00 

CM 

o 

d 

05 

o 

d 

CD 

o 

d 

X — 

o 

d 

o 

o 

o 

o 

d 

05 

o 

o 

d 

1 

CO 

LO 

05 

CM 

CM 

00 

d 

CM 

Is- 

O; 

CO 

Is- 

d 

00 

r- 

00 

xt 

d 

oo 

LO 

LO 

CD 

CD 

o 

05 

CO 

CD 

CO 

CM 

q 

CO 

1 

L— 

0 

0 

E 

TO  (D 
0  2 
CL  0 
<  > 

o 

05 

X — 

o 

00 

o 
1^ 
x — 

o 

CD 

X- 

o 

LO 

o 

X — 

o 

oo 

o 

CM 

x— 

o 

o 

o 

o 

LO 

o 

i— 

0 
-»— > 

0 

i  ® 

E  = 
0  0 
CL  > 

o 

o 

CM 

o 

00 

o 

o 

LO 

o 

CD 

o 

r- 

o 

00 

o 

05 

o 

o 

o 

LO 

o 

o 

CM 

O 

LO 

CM 

Q 

< 

ULI 

< 

Is- 

LO 

LO 

00 

CM 

00 

Is- 

LO 

00 

CD 

LO 

00 

o 

LO 

CM 

05 

Is- 

00 

05 

CD 

LO 

05 

CD 

Is- 

CD 

Is- 

CM 

CD 

CD 

X — 

Is- 

CD 

r- 

o 

o 

Is- 

CD 

CD 

TT 

c5 

05 

CD 

CD 

CD 

o 

d 

05 

Is- 

CD 

CD 

o 

Is- 

CD 

CD 

CD 

05 

00 

CM 

CD 

CD 

CD 

CD 

o 

d 

00 

LO 

CD 

CD 

00 

LO 

CO 

o 

CD 

CD 

o 

o 

5 

< 

LU 

< 

Is- 

LO 

LO 

00 

CM 

00 

Is- 

LO 

CM 

05 

CD 

CD 

o- 

05 

CD 

CD 

LO 

CD 

LO 

X — 

CM 

CM 

o 

LO 

CD 

1^ 

00 

00 

00 

Is- 

CD 

s 

d 

o 

LO 

3 

LO 

05 

05 

Is- 

00 

CD 

00 

oo 

LO 

X— 

CO 

CM 

05 

CM 

CM 

CD 

LO 

CM 

d 

o 

CD 

CD 

CM 

d 

oo 

05 

00 

CD 

05 

Is- 

1^ 

05 

Is- 

co 

CD 

o 

LO 

Q 

< 

LU 

< 

00 

CD 

CD 

r- 

LO 

CD 

CM 

d 

LO 

LO 

CM 

CD 

Is- 

CO 

Is- 

Is- 

5 

Is- 

00 

00 

Is- 

co 

s 

X - 

d 

00 

o 

T~ 

CD 

o 

CM 

LO 

LO 

05 

o 

CD 

CO 

o 

d 

00 

o 

CD 

05 

LO 

1^ 

CD 

o 

CD 

LO 

d 

CD 

CD 

o 

CD 

00 

00 

d 

CM 

LO 

o 

CD 

CO 

d 

oo 

xt 

o 

CD 

05 

X— 

CO 

o 

CD 

D 

< 

LU 

LO 

h- 

o 

Is— 

r- 

CD 

q 

LO 

LO 

CM 

Is- 

o 

Is-; 

LO 

o 

Is- 

LO 

q 

00 

CD 

o 

o 

r- 

X — 

LO 

05 

r— 

05 

CD 

o 

d 

o 

CD 

05 

CD 

CO 

05 

d 

3 

CD 

CM 

Is-' 

CO 

05 

CD 

Xtf 

rj- 

CO 

oo 

CM 

05 

CD 

00 

00 

d 

00 

o 

05 

CD 

CM 

Is- 

r- 

00 

00 

00 

CD 

26 


MinE  0.009  MinE50  0.QQ03  MinEIOQ  0.007  MinE250  0.Q07 
MaxE  0.125  MaxE50  Q.093Q  MaxEIOQ  0.157  MaxE250  0.109 


Output  File  -  Sensitivity  to  Uncertainty 
Flow  Frequency  Function 


Ji i  ^  rv* 


O) 

c 

<D  > 

__  CD 

-O  -D 

O  “O 
o  ^ 
cd  to 

-  CD 
CD  U) 

E>  CO 
_TO  — 

cr  § 
P  c 


CD 

£Z 

CD 


O 

O 

0 


>* 

■o 

3 

co 

.  CD 
>  w 
CD 

^  O 


CO  £ 

2  o  £ 

c  U-  CO 

3 

O' 

15 

.y  z  c/5 

a  -i  -J 

>  u_  Li- 
H  LL  C/5 


CD 

c  S 

w  P 

„  Q_ 
<D  O 
CD  £ 

S  c/5 

CO 

k  ® 

CD 

CO  t+T 


0 

CL 

O 


CO  _ 
CO  CO 
CO  Li! 
Q  LL. 


0 

Q_ 

O 

CO 


0 

E 

TO 


C 

3 

a: 


o 

LO 

N~ 

o 

CD 

00 

CO 

CO 

LO 

h- 

CM 

O 

CD 

LO 

LO 

M- 

M- 

CO 

CO 

CNJ 

00 

CM 

X — 

o 

o 

o 

O 

p 

P 

o 

o 

o 

CO 

LU 

d 

d 

d 

d 

d 

d 

d 

d 

o 

o 

_d 

El 

d 

o 

m- 

CD 

CD 

CO 

CO 

CO 

h- 

o 

o 

o 

o 

LO 

CM 

LU 

o 

LO 

o 

m- 

LO 

X — 

M* 

CM 

XT“ 

o 

o 

00 

CM 

o 

CD 

LO 

CM 

CM 

X — 

X — 

X — 

X — 

X — 

t — 

o 

LU 

LO 

O 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

X 

CM 

O 

p 

O 

o 

o 

o 

o 

p 

o 

o 

o 

■ 

c 

TO 

CO 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

■ 

o 

m- 

CO 

CD 

CD 

o 

CD 

CD 

CM 

CM 

N- 

r- 

■ 

N- 

M’ 

o 

CO 

00 

T— 

CO 

r^- 

LO 

M" 

M* 

M" 

00 

CM 

CM 

CD 

X — 

CO 

T- 

o 

o 

p 

o 

O 

O 

p 

o 

O 

CO 

uu 

d 

o 

d 

o 

d 

d 

d 

d 

d 

d 

d 

■ 

o 

o 

CD 

o 

LO 

CD 

LO 

CM 

X"~ 

CD 

CD 

00 

CD 

o 

o 

o 

o 

o 

r^ 

CM 

X — 

X — 

X — 

o 

o 

o 

O 

LU 

UJ 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

X 

o 

o 

o 

o 

o 

p 

p 

p 

o 

p 

O 

TO 

CO 

o 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

CO 

00 

CD 

CO 

N- 

CO 

LO 

o 

CD 

CM 

CM 

CO 

o 

CO 

CD 

o 

00 

CD 

LO 

M- 

M- 

CO 

CO 

CO 

CO 

CD 

LO 

CO 

T— 

T— 

o 

o 

o 

o 

O 

o 

o 

o 

m 

o 

CO 

1X1 

o 

d 

d 

d 

d 

o 

d 

d 

o 

d 

d 

m 

o 

d 

r^- 

CD 

CD 

r^ 

LO 

LO 

CO 

O 

o 

00 

X — 

o 

LO 

CM 

00 

r- 

h*» 

CD 

LO 

LU 

LO 

CD 

CO 

CM 

X — 

X — 

X- 

o 

o 

o 

o 

o 

LU 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

X 

LO 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

TO 

CO 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

00 

CD 

CD 

h- 

T- 

CD 

00 

o 

o 

00 

CD 

N- 

T— 

00 

h- 

LO 

M- 

M" 

M" 

CO 

CO 

■ 

CO 

CD 

CO 

t— 

T— 

p 

p 

p 

O 

p 

O 

p 

o 

1 1 

o 

CO 

Lit 

d 

o 

d 

d 

d 

d 

o 

d 

d 

d 

d 

o 

d 

h- 

h- 

T— 

CO 

CD 

LO 

CO 

CM 

o 

00 

| 

LU 

O) 

M* 

00 

CM 

T— 

X — 

X — 

X — 

T“ 

X — 

o 

LU 

o 

O 

o 

o 

o 

o 

o 

o 

o 

O 

o 

X 

o 

O 

o 

o 

o 

o 

o 

o 

o 

O 

o 

■tv 

_co 

o 

d 

d 

d 

o 

o 

d 

d 

d 

d 

d 

■ 

S— 

0 

0 

E 

CO 

0 

TO 

_3 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

CL 

TO 

CD 

00 

N- 

CD 

LO 

M- 

CO 

CM 

x— 

o 

o 

< 

> 

T” 

■*“ 

x — 

LO 

o 

0 

0 

E 

0 

CO 

13 

l_ 

o 

o 

o 

TO 

"to 

o 

o 

o 

o 

o 

O 

o 

o 

O 

o 

LO 

o 

CL 

_> 

X — 

CM 

CO 

M" 

LO 

CD 

N- 

00 

CD 

x- 

x— 

CM 

O 

LO 

CM 

a 

< 

CD 

CO 

CD 

o 

LO 

X~ 

CD 

CD 

LO 

CD 

LO 

Hi 

LO 

00 

CM 

CM 

is 

o 

o 

O) 

CD 

CD 

p 

co 

CD 

LO 

LO 

LO 

EE 

LO 

LO 

M1 

M- 

M" 

o 

o 

T — 

D 

< 

00 

CD 

o 

N. 

o 

M- 

X — 

00 

N- 

LO 

O 

h* 

LU 

00 

O 

00 

p 

p 

p 

LO 

M; 

M; 

M; 

M* 

CO 

<_ 

LO 

LO 

M- 

M- 

M- 

M1 

M- 

M- 

o 

LO 

D 

< 

h- 

M- 

CM 

CM 

CD 

x — 

00 

CD 

LO 

M* 

O 

h- 

LU 

q 

CO 

T“ 

p 

CD 

CD 

00 

00 

CO 

CO 

00 

N- 

<_ 

LO 

M- 

M- 

M- 

CO 

CO 

CO 

CO 

CO 

CO 

co 

CO 

Q 

m- 

M- 

CD 

h- 

o 

CO 

CO 

CD 

CO 

CD 

< 

T“ 

T— 

00 

CD 

LO 

LO 

M" 

M; 

M; 

CO 

CO 

CM 

LU 

CD 

LO 

LO 

LO 

LO 

LO 

LO 

LO 

LO 

LO 

UD 

127 


Appendix  A 

Uncertainty  Analysis  Results  -  EAD 

D2575  =  Difference  between  the  25%  and  75%  exceedance  values  in  EAD  distribution 
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Uncertainty  Analysis  Results  -  EAD  With  50-yr  Levee 

D2575  =  Difference  between  the  25%  and  75%  exceedance  values  in  EAD  distribution 
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S  EAD95  EAD75  EAD25  EAD5  D2575  S  IEAD95  IEAD75  IEAD25  IEAD5  ID2575 

Base  0.00481  0.00712  0.01160  0.01611  0.00448  Base  0.00477  a00709  1 0.01 157  0.01612  0.00448 

5%  0.00479  0.00708  0.01162  0.01617  0.00454  5%  0.00474  0.00706  0.01157  0.01619  0.00451 

10%  0.00468  0.00700  0.01167  0.01643  0.00467  10%  0.00464  0.00698  0.01166  0.01639  0.00468 

20%  1 0.00426  10.00667  1 0.01 175  10.01717  1 0.00508  1  1 20%  1 0.00426  1 0.00668  10.01178  1 0.01 722  1 0.00510 


Uncertainty  Analysis  Results  -  EAD  With  250-Yr  Levee 

D2575  =  Difference  between  the  25%  and  75%  exceedance  values  in  EAD  distribution 
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S  EAD95  EAD75  EAD25  EAD5  D2575  S  IEAD95  IEAD75  IEAD25  IEAD5  ID2575 

Base  0.00177  0.00273  0.00491  0.00734  0.00218  Base  0.00177  0.00275  0.00493  0.00740  0.00218 

5% _ 0.00175  0.00272  0.00491  0.00738  0.00219  5% _  0,00177  0.00274  0.00494  0.00745  0.00220 

10%  0.00173  0.00269  0.00493  0.00745  0,00224  10%  0,00174  0,00270  0,00495  0.00750  0.00225 

20%  1 0.001 59  1 0.00259  1 0.00499  1 0,00778  1 0.00240  120%  Iq.00160  Iq.00260  10.00501  lp.00784  Iq.00241 


